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FOREWORD 


The  work  reported  herein  was  sponsored  in  part  by  the 
United  States  Army  Research  Office,  Durham,  under  Contract 
No.  DA-31-124-ARO--D-246 .  The  study  presented  herein  was 
conducted  by  Mr.  Ronald  A.  Cudnick  in  fulfilling  the  require¬ 
ments  for  a  thesis  in  his  Master  of  Science  program  at  The 
Ohio  State  University. 


ABSTRACT 


A  study  was  made  of  charged  particle  notion  in  a  frec-vortex 
flow  field  to  determine  the  parameters  affecting  particle  motion 
and  to  determine  the  extent  to  which  applied  electric  fields  can 
influence  the  particle  motion. 

Four  different  cases  were  investigated.  These  included  first, 
the  analysis  of  the  motion  of  an  uncharged  particle  in  a  free-vortex; 
second,  analysis  of  the  motion  of  a  charged  particle  in  a  viscous 
medium  under  the  influence  of  an  applied  electrostatic  field; 
third,  analysis  of  charged  particle  motion  in  a  free-vortex  under 
the  influence  of  the  applied  electrostatic  field,  assuming  the 
particle  is  first  positively  and  then  negatively  charged;  and  fourth, 
analysis  of  the  motion  of  two  charged  particles  in  a  free-vortex, 
accounting  for  field  effects  due  to  particle  charge. 

The  results  indicate  that  in  all  cases,  as  the  particle 
diameter  and  density,  and  free-vortex  flow  velocity  are  increased, 
radial  particle  migration  is  increased.  Radial  particle  migration 
is  enhanced  if  the  particle  is  assumed  to  be  positively  charged  and 
under  the  influence  of  the  electrostatic  field  due  to  a  positively 
charged  line  source.  If  the  particle  is  assumed  to  be  negatively 
charged,  the  electrostatic  force  tends  to  balance  the  centrifugal 
force  and  allow  the  particle  to  maintain  an  equilibrium  radius. 
Finally,  the  field  effects  due  to  the  particle  charges  themselves 
are  insufficient  to  affect  particle  motion  in  the  two  particle  system. 
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INTRODUCTION 


The  motion  of  charged  and  uncharged  particles  in  a  free- 
vortex  flow  field  is  of  interest  in  a  wide  variety  of  areas.  These 
areas  range  from  particle  separation  and  classification  by  centrifugal 
and  electrostatic  action,  to  possibly  improved  operation  of  electro¬ 
static  precipitators  by  addition  of  some  device  to  impart  swirling 
motion  to  the  exhaust  gases,  to  the  extremely  interesting  area  of 
the  role  of  charged  particles  in  the  initiation  and  maintenance  of 
two  of  nature's  vortex  systems  --  the  tornado  and  the  dust  devil. 

When  this  thesis  was  first  undertaken,  it  was  with  the  goal  of 
gaining  some  insight  into  determining  the  role  played  by  the 
electrical  activity  associated  with  tornadoes  and  dust  devils  in 
the  initiation  and  maintenance  of  these  vortex  systems.  While  it 
was  soon  realized  that  this  was  too  complex  a  task  to  be  undertaken 
for  a  thesis,  the  problem  was  still  a  fascinating  one,  and  as  a 
result  a  literature  survey  was  made  to  determine  the  current  theories 
on  the  significance  of  the  electrical  activity  associated  with 
tornadoes  and  dust  devils.  The  information  gained  in  the  literature 
survey  has  been  summarized  and  is  discussed  in  the  next  section. 

This  survey  of  the  literature  also  indicated  that  while  work 
had  been  done  in  determining  the  motion  of  uncharged  particles  in  a 
swirling  flow  by  Hirsthkron  and  Ehrich  (Ref.  1),  Kriebel  (Ref.  2), 
Lapple  and  Shepherd  (Ref.  3)  and  Uematu,  et  al.,  (Ref.  4),  no  woi.1 
had  been  done  in  determining  the  motion  of  charged  particles  in  a 
swirling  flow  under  the  influence  of  an  applied  electrostatic  field. 


It  was  then  decided  that  the  objective  of  this  thesis  would  be  to 
study  the  motion  of  a  charged  particle  in  a  free-vortex  under  the 
influence  of  an  applied  electric  field  and  to  determine  the  parameters 
influencing  the  particle  motion. 

The  significance  of  electrification  on  the  dynamics  of  a 
particulate  system  has  been  extensively  studied  by  Soo  (Ref.  5,  6, 

7,  and  8).  His  work  has  shown  that  electrification  of  the  solid 
particles  in  a  gas-solid  suspension  affects  the  concentration 
distribution  of  the  solid  particles  and  would  be  expected  to  affect, 
for  example,  the  friction,  heat  transfer,  and  settling  of  solid 
particles.  However,  his  work  has  been  primarily  concerned  with  one¬ 
dimensional  flow  through  pipes  rather  than  swirling  flow. 
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BACKGROUND  OF  THE  ROLE  OF  ELECTRICAL  ACTIVITY 
ASSOCIATED  WITH  VORTEX  SYSTEMS,  INCLUDING 
TORNADOES  AND  DUST  DEVILS 

Determining  the  role  o£  electrical  activity  in  vortex  systems, 
in  particular  tornadoes  and  dust  devils,  is  a  complex  problem 
which  has  puzzled  man  for  many  centuries.  The  idea  that  electrical 
activity  may  be  a  cause  rather  than  a  result  of  tornadoes  is 
certainly  not  new,  as  evidenced  by  the  works  of  Lucretius  (60  B.C.) 
and  Francis  Bacon  (1622)  cited  by  Vonnegut  (Ref.  9)  in  his  dis¬ 
cussion  of  the  electrical  theory  of  tornadoes.  The  various  accounts 
of  the  electrical  phenomena  associated  with  tornadoes,  dust  devils, 
and  confined  vortex  systems  gives  much  credence  to  the  belief  that 
the  electrical  activity  and  the  vortex  motion  are  more  than 
fortuitously  interrelated.  These  accounts  also  serve  as  the  basis 
for  theories  proposed  to  explain  the  role  of  electrical  activity 
in  vortex  systems. 

The  accounts  of  electrical  activity  associated  with  vortex 
systems  to  be  presented  here  are  certainly  not  all  that  are  available 
it  is  hoped,  however,  that  these  accounts  will  give  a  representative 
sampling  of  the  type  of  electrical  behavior  that  has  been  observed 
in  tornadoes,  dust  devils,  and  confined  vortex  systems. 

Freier  (Ref.  10)  obtained  measurements  of  the  earth's  electric 
field  in  the  Minneapolis- St .  Paul  area  while  tornadoes  were  within 
a  seventy-mile  radius  of  the  recording  station.  The  measurements 
he  obtained  show  a  high  frequency  oscillation  of  the  earth's  electric 
field  beginning  nearly  at  the  time  of  the  first  reported  tornado 
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and  continuing  marly  an  hour  after  the  last  tornadoes  struck.  The 
significant  points  noted  from  these  measurements  were  first,  that  the 
high  frequency  oscillations  of  the  earth's  electric  field  present  had 
not  been  observed  at  the  recording  station  during  any  other  thunder¬ 
storm  activity  during  the  previous  year  and  second,  the  measurements 
were  obtained  when  nowhere  within  the  visible  horizon  was  there  an 
indication  of  a  storm.  Freier  states  that  while  it  is  difficult  to 
draw  any  conclusions  from  the  limited  data,  it  is  significant  that 
the  high  frequency  oscillations  of  the  earth's  electric  field  are 
observed  only  during  the  tornado  activity. 

Vonnegut  and  Weyer  (Ref.  11)  confirmed  the  existence  of  unusual 
electrical  discharges  in  a  tornado  at  Toledo,  Ohio  on  April  11,  1965. 
They  present  a  photograph  showing  two  illuminated  vertical  pillars 
indicating  the  presence  of  some  unusual  luminous  phenomena.  From 

t 

eyewitness  reports  and  the  paths  of  destruction,  it  appeared  that 
there  were  two  tornadoes  following  essentially  the  same  path.  Over 
a  dozen  eyewitness  accounts  were  obtained,  many  of  which  indicated 
various  kinds  of  luminosity.  Many  of  the  eyewitnesses  reported  a 
white  or  i  luish-wh  i L e  light,  definitely  not  lightning,  but  sufficiently 
bright  to  illuminate  the  surroundings  and  give  the  appearance  of  day¬ 
light.  There  were  no  observations  of  any  significant  lightning 
strokes.  Several  of  the  eyewitnesses  observed  bright  balls  of  light, 
either  orange  or  multicolored.  These  observations  suggest  the 
presence  of  ball  lightning,  a  phenomenon  of  which  very  little  is 


known . 


In  discussing  the  two  columns  of  light  appearing  in  the 
photograph,  the  authors  suggest  that  the  columns  of  light  may  be 
due  to  some  kind  of  luminous  electrical  discharge  or  they  may  be 
tornado  funnels  illuminated  by  lightning  or  some  other  type  of  dis¬ 
charge  within  the  vortex.  While  there  are  no  other  photographs  of 
a  luminous  tornado  funnel,  there  have  been  several  other  observations 
of  an  apparently  similar  phenomenon. 

Brook  (Ref.  12)  reported  a  measurement  of  the  perturbations  in 
the  magnetic  field  and  the  earth  current  taken  in  the  vicinity  of 
a  tornado.  The  data  were  obtained  by  Boucher  near  Tulsa,  Oklahoma 
on  May  27,  1962.  The  significant  feature  of  the  measurements 
made  is  a  large  rapid  deflection  of  the  horizontal  earth  current. 
Coincident  with  this  deflection  was  a  large  sudden  shift  in  the 
horizontal  component  of  the  magnetic  field,  with  smaller  deflections 
in  the  declination  and  vertical  component.  The  coincident,  step¬ 
like  deflections  of  all  the  parameters  measured  raise  the  possibility 
that  a  common  source  existed  for  all  the  disturbances.  Brook 
believes  that  the  source  of  the  disturbance  was  the  tornado  that 
was  seen  touching  down  at  the  precise  time  at  which  the  deflections 
occurred . 

Vonnegut  and  Moore  (Ref.  13)  discuss  electrical  field  measure¬ 
ments  made  by  Gunn  (Ref.  14)  near  a  tornado  on  May  25,  1955  at 
Blackwell,  Oklahoma  and  Udall,  Kansas.  In  his  work  Gunn  concluded 
that  the  measurements  give  little  evidence  that  first,  electrical 
effects  near  the  tornado  differ  basically  from  ordinary  thunderstorm 


electrification  and  second,  that  the  tornado  is  primarily  an 
electrical  phenomenon.  Vonnegut  and  Moore  postulate  that  although 
electrical  effects  are  more  severe  in  tornadoes  than  in  ordinary 
thunderstorms,  they  are  basically  not  different.  Thus,  they  agree 
with  Gunn's  opinion  that  the  electric  field  measurements  made  near 
the  tornado  are  what  would  be  expected  in  ordinary  thunderstorms  if 
an  order  of  magnitude  increase  in  the  turbulent  velocity  occurred. 
Eyewitnesses  to  the  tornado  in  Blackwell  observed  evidence  of  intense 
electrical  activity  similar  to  that  reported  by  Vonnegut  and  Weyer 
in  the  Toledo  tornado.  According  to  the  authors,  the  luminous  dis¬ 
charges  in  a  tornado  suggest  that  large  current  densities  and 
electric  fields  exist  inside  the  funnel  and  on  the  ground  directly 
beneath  it. 

Silberg  (Ref.  15)  reports  that  there  have  been  repeated  ob¬ 
servations  of  a  tornado  burning,  scorching,  or  dehydrating  trees 
and  vegetation  in  its  path.  In  the  past  these  heating  effects  were 
explained  with  the  aid  of  a  "hot  wind"  or  a  corona  discharge.  Silberg 
postulates  that  an  intense  ring  current  existing  in  the  tornado  cloud 
could  be  responsible  for  these  heating  effects  and  ocher  phenomena 
associated  with  a  tornado.  The  theory  Silberg  presents  for  this  will 
be  discussed  in  greater  detail  in  another  section  of  this  summary. 

Silberg  (Ref.  16)  also  presents  data  obtained  from  passive 
electrical  measurements  from  three  tornadoes  in  Oklahoma  during 
the  spring  of  1963.  The  tornado  measurements,  when  compared  with 
those  of  a  local  thunderstorm,  show  that  the  electrical  discharges 
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within  the  tornadoes  are  much  more  intense  than  in  the  local 
thunderstorm,  particularly  in  the  frequency  of  occurrence. 

Several  observers  have  obtained  measurements  of  the  electric 
field  of  a  dust  devil.  Freier  (Ref.  17)  obtained  a  record  of  the 
electric  field  of  a  dust  devil  in  the  Sahara  Desert.  The  record 
shows  that  the  electric  field  decreased  from  its  average  value 
of  about  150  volts  per  meter  before  and  after  passage  of  the  dust 
devil  to  -75  volts  per  meter,  increased  to  approximately  200  volts 
per  meter  and  then  decreased  again  to  -450  volts  per  meter  before 
returning  to  its  average  value.  Freier  compares  the  experimental 
data  to  that  predicted  by  a  simple  dipole  model  with  negative  charge 
above . 

Crozier  (Ref.  18)  recorded  the  electric  field  for  a  fairly 
large  dust  devil  passing  within  1500  feet  of  the  electric  field 
recording  instrument.  The  potential  gradient  record  indicated  that 
the  gradient  was  lowered  for  approximately  10  minutes  and  the  lowest 
value  was  approximately  100  volts  per  meter  below  the  average  value 
before  and  after  passage  of  the  dust  devil.  Crozier  presents  a 
simple  model  to  explain  the  potential  gradient  record  for  the  dust 
devil  and  also  suggests  mechanisms  for  obtaining  the  space  charge 
densities  necessary  in  his  simple  electrostatic  model.  Crozier  attempts 
to  compare  the  electrical  properties  of  the  dust  devil  he  observed  to 
those  of  the  Sahara  dust  devil  described  by  Freier  (Ref.  17).  Crozier 
suggests  that  the  W-shape  of  the  Sahara  dust  devil  gradient  record 
could  be  indicative  of  a  rapid  increase  in  electrical  activity  within 
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the  dust  devil.  This  model  will  be  discussed  in  detail  in  a  following 
section , 

Bradley  and  Semonin  (Ref.  19)  conducted  airborne  electrical 
measurements  in  dust  whirls  in  Illinois.  Measurements  were  made  of 
the  vertical  potential  gradient  and  space  charge  concentration  in  the 
tops  of  three  large  dust  whirls.  In  all  cases  the  measurements  showed 
the  presence  of  a  negative  space  charge  with  an  increase  in  the 
potential  gradient  indicating  an  excess  of  negative  charge  below. 

Silberg  and  Goshgarian  (Ref.  20)  described  several  experiments 
which  demonstrate  that  rotary  and  vortical  motion  can  be  produced  by 
initially  electrostatic  fields.  A  circular  array  of  16  electrodes 
was  installed  in  a  large  lucite  tube.  Two  adjacent  electrodes, 
sufficiently  separated  to  prevent  complete  breakdown,  were  connected 
to  a  90  kilovolt  power  supply.  A  small  smoke  pot  was  placed  on  the 
tube  axis  three  feet  down  from  the  electrode  array  and  activated. 

When  the  power  supply  was  turned  on,  votary  fluid  motion  was  started 
and  continued  as  long  as  the  power  was  on. 

Later  work  by  Silberg,  Goshgarian,  and  Johnson  (Ref.  21) 
attempted  to  measure  the  electrostatically  produced  vortical  motion. 
The  equipment  used  was  the  same  with  the  exception  of  a  larger  power 
supply  ( 1 20  kilovolts).  Air,  bromine  vapor,  and  smoke  mixed  with 
air  were  the  fluids  used.  No  vortical  motion  was  observed  at  the  power 
supply  voltages  used  for  air  and  for  bromine  vapor.  Using  smoke  mixed 
with  air,  at  30  kv  considerable  turbulence  was  noted  around  the 
electrodes.  Above  50  kv ,  the  smoke  particles  at  the  top  and  bottom 
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assumed  a  regular  rotary  motion  and  in  a  short  time  the  entire 
smoke  column  began  rotating,  increasing  uniformly  in  velocity  until 
a  well  defined  vortex  was  formed.  An  attempt  was  made  to  measure 
the  fluid  motion  using  a  Pitot-static  tube  with  a  minimum  detectible 
velocity  calibrated  to  be  approximately  10  feet  per  second.  When  a 
traverse  of  the  Pitot  tube  was  made  through  the  vortex  about  4  inches 
above  the  mouth  of  the  smoke  pot  where  the  smoke  velocity  appeared  to 
be  greatest,  no  velocity  was  measured.  It  appeared  that  if  there  was 
motion  of  the  air  caused  either  by  the  electrostatic  field  or  by  drag 
induced  by  the  accelerated  smoke  particles,  then  this  motion  was  below 
the  minimum  detectible  threshold  of  the  Pitot  tube.  Several  con¬ 
clusions  were  reached  by  the  authors:  first,  the  observed  vortical 
motion  existed  primarily  in  the  dielectric  smoke;  second,  the  vortical 
motion  in  the  smoke  was  high,  although  immediately  exterior  to  it 
no  measurable  air  velocity  was  detectible;  third,  from  the  high 
potentials  required  to  initiate  vorticity,  it  was  concluded  that  the 
dielectric  particles  to  which  the  field  initially  couples  were 
probably  charged;  and  fourth,  since  the  vortical  motion  did  not 
extend  to  the  boundary,  no  boundary  effects  were  noted. 

Lavan  and  Fejer  (Ref.  22)  observed  a  luminescent  region  near  the 
axis  of  a  cylindrical  duct  confining  a  supersonic  swirling  air  flow 
at  ambient  temperature.  This  luminescence  was  identified  as  a  glow 
discharge  induced  by  a  strong  internal  electric  field  produced  by 
condensed  water  droplets  impinging  on  the  duct  wall.  It  was  found 
that  the  luminescence  increased  uniformly  with  swirl  intensity  and 
with  the  specific  humidity  of  the  ambient  air,  while  it  decreased 
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when  the  axial  Mach  number  was  increased  or  decreased.  When  the 
dielectric  duct  was  replaced  by  one  made  of  a  conducting  material,  no 
luminescence  was  observed,  suggesting  that  electrostatic  charges  on 
the  surface  of  the  dielectric  wall  were  a  likely  cause  of  the 
phenomenon. 

The  various  accounts  presented  here  strongly  suggest  that  an 
interrelationship  does  exist  between  electrical  activity  or  phenomena 
and  vortex  systems.  Several  theories  have  been  proposed  to  explain 
the  role  of  electrical  activity  in  the  initiation  and  continuation 
of  tornadoes.  The  two  main  theories  here  are  the  electrical  heating 
theory  and  the  electromagnetic  or  electrodynamic  initiation  theory. 

In  these  theories  the  electrical  activity  is  the  cause  of  tornadoes 
rather  than  a  result  of  them.  Other  theories  suggest  that  the 
electrical  activity  is  the  result  of  vortex  motion.  The  various 
theories  will  be  presented  and  briefly  discussed  here. 

Vonnegut  (Ref.  9)  suggests  that  there  is  sufficient  electrical 
energy  in  an  intense  thunderstorm  to  power  a  tornado  and  that  the 
electrification  could  cause  very  intense  winds  by  electrically  heating 
air  or  by  accelerating  charged  air  in  an  electric  field.  Vonnegut 
notes  that  to  obtain  winds  of  the  velocity  that  appear  to  exist  in  a 
tornado  there  must  be  temperature  differences  between  air  masses  in 
the  atmosphere  probably  of  the  order  of  several  hundred  degrees 
centigrade.  He  concludes  it  is  doubtful  that  ordinary  atmospheric 
processes  could  produce  the  large  temperature  differences  required 
for  winds  of  tornado  speed.  He  notes  that  the  quantities  of  heat 
sufficient  to  raise  large  volumes  of  air  to  a  high  temperature  are 


released  by  volcanic  activity  or  large  fires  and  that  tornadoes  or 

fire  whirls  have  been  observed  under  these  conditions.  His  electrical 

theory  of  tornadoes  uses  as  its  bases  observations  of  tornadoes  which 

indicate  the  presence  of  sufficient  electrical  energy  to  account  for 

their  intensity,  and  the  existence  of  mechanisms  by  which  this  electrical 

activity  could  produce  the  intense  winds.  The  accounts  of  electrical 

activity  in  tornadoes  presented  previously  will  serve  to  verify  its 
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presence  and  Vonnegut  estimates  that  roughly  10  kw  of  electrical 
power  is  available  in  a  tornado  producing  storm.  It  thus  appears 
that  sufficient  electrical  energy  probably  exists  in  a  violent 
thunderstorm  to  cause  a  tornado,  if  the  energy  is  properly  applied. 
Vonnegut  next  attempts  to  establish  that  this  electrical  energy  can 
be  transformed  into  the  kinetic  energy  of  rapidly  moving  air.  Two 
mechanisms  by  which  this  can  happen  are  presented.  The  first  is 
that  electrical  energy  could  heat  a  volume  of  air  to  such  a  high 
temperature  that  very  intense  convection  would  take  place.  The  second 
is  that  highly  charged  air  could  be  accelerated  to  a  high  velocity  in 
a  strong  electric  field.  Both  processes  are  dependent  on  the  trans¬ 
port  of  electric  charge  through  a  potential  difference.  Vonnegut 
considers  three  charge  transport  mechanisms  that  might  operate  in  a 
storm. 

The  first  mechanism  is  the  lightning  stroke,  in  which  a  narrow 
path  through  the  air  suddenly  becomes  highly  ionized  and  conducting 
with  the  charge  flowing  in  a  fraction  of  a  second.  This  is  an 
intermittent  process,  however,  in  which  the  instantaneous  current 
may  reach  many  kiloamperes,  although  the  average  current  is  much  less. 
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The  second  mechanism  is  a  continuous  arc  or  glow  discharge  in 
which  current  may  flow  and  heat  the  air.  Here  the  current  flow  is 
steady  and  does  not  reach  the  high  instantaneous  values  of  the 
lightning  discharge. 

The  third  mechanism  is  by  the  movement  of  charged  particles  such 
as  ions  or  electrified  cloud  particles  in  an  electric  field.  These 
conduction  processes  probably  account  for  the  largest  part  of  the 
charge  transport  from  the  earth  to  the  clouds  in  thunderstorms. 

In  all  three  processes  a  large  portion  of  the  electrical  energy 
is  released  as  heat.  Vonnegut  suggests  that  each  of  these  three 
mechanisms  can,  under  the  right  circumstances,  cause  intense  winds 
by  convection  resulting  from  the  hot  air  produced.  The  first  two 
mechanisms  are  dependent  on  the  efficiency  of  the  conversion  of  energy 
released  by  lightning  in  the  form  of  heat  into  kinetic  energy  of  the 
wind.  A  simple  calculation  shows  that  an  air  mass  could  be  con¬ 
tinuously  heated  by  more  than  100 °C  under  these  conditions. 

In  the  third  mechanism,  the  flow  of  charge  by  conduction, 
transformation  of  electrical  energy  into  kinetic  energy  can  take  place 
directly  if  the  charge  transport  occurs  by  the  movements  of  ions  all 
having  predominantly  the  same  sign.  The  movement  of  these  ions  under 
the  influence  of  an  electric  field  will  transfer  momentum  directly  to 
the  air  and  cause  it  to  accelerate.  The  flow  of  current  will  cause 
it  to  move  the  air  and  heat  it  under  the  influence  of  electrical  forces. 
The  ratio  of  the  fraction  of  the  electrical  energy  transformed  into 
hear  to  that  transformed  into  mechanical  work  will  be  equal  to  the 
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ratio  between  the  velocity  of  the  ions  relative  to  the  air  and  the 
velocity  of  the  air  itself.  The  importance  of  the  electrical 
acceleration  of  the  air  will  be  greatest  when  the  wind  velocity  is 
high  and  the  mobility  of  the  charge  carriers  is  low.  The  mobility 
will  be  smallest  when  the  charge  is  carried  predominantly  on  dust  or 
cloud  particles.  Vonnegut  presents  a  simple  model  showing  that  the 
electrical  force  acting  on  a  charged  water  droplet  in  a  typical 
thunderstorm  electrical  field  is  equivalent  to  that  force  which  would 
result  if  the  air  were  heated  20 °C  above  its  environment.  Thus  it 
would  appear  that  the  first  two  charging  mechanisms  would  be  primarily 
responsible  for  any  vortex  motion  due  to  electrical  heating.  The 
third  mechanism  would  contribute  to  the  electrical  heating  effect, 
but  its  more  important  role  may  be  that  of  accelerating  the  air  due 
to  the  movement  of  ions  under  the  influence  of  an  electrical  field. 

Wilkins  (Ref.  23)  conducted  laboratory  model  vortex  experiments 
and  theoretical  investigations  to  determine  whether  the  unusual 
electrical  phenomena  sometimes  accompanying  tornadoes  might  play  an 
important  role  in  the  initiation  or  continuation  of  a  tornado.  Wilkins 
concluded  that  first,  heat  released  by  lightning  may  serve  to  initiate 
updrafts,  but  its  effectiveness  in  tornado  formation  is  difficult  to 
evaluate  due  to  the  presence  of  the  latent  heat  available  in  a  con- 
vectively  unstable  atmosphere.  Second,  the  electrodynamic  accelerations 
of  ions  appear  to  be  unimportant  in  either  the  formative  or  mature 
stages  of  a  tornado. 

Colgate  (Ref.  2A)  postulates  that  the  tornado  vortex  is  driven 
by  a  line  sink  of  electrically  heated  air  extending  at  least 
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5  kilometers  high.  Ilis  argument  for  rn  electrical  energy  mechanism 
is  that  purely  hydrodynamic  motion  cannot  adequately  account  for  the 
very  high  tornado  wind  speeds  derivable  from  a  maximum  available 
adiabatic  pressure  difference  of  a  tenth  of  an  atmosphere  from  the 
ground  to  the  base  of  the  tropopause. 

Theories  proposed  by  Vonnegut  (Ref.  9),  Rathbun  (Ref.  25),  and 
Lewellen  (Ref.  26)  suggested  that  a  tornado  or  a  vortex  could  be 
electrodynamically  or  clectromagnetically  driven.  The  theory  proposed 
by  Vonnegut  on  the  acceleration  of  air  by  the  motions  of  ions  under  the 
influence  of  an  electric  field  has  already  been  discussed  and 
Wilkins'  (Ref.  23)  objections  to  it  have  been  stated.  Rathbun  (Ref.  25) 
proposed  an  electromagnetic  basis  for  the  iniLiaLion  of  a  tornado.  lie 
suggests  that  the  process  leading  to  the  formation  of  a  tornado  is 
the  following: 

(1)  Electric  charge  accumulation  in  a  cloud  increases 
until  the  atmosphere  is  almost  stressed  to  the*  point 
of  breakdown. 

(2)  A  leader  stroke  of  lightning  develops,  but  because 
of  the  energy  available,  it  is  not  followed  by  a 
full  lightning  stroke.  This  threatening  stroke  is 
believed  to  exist  lor  several  seconds. 

(3)  A  considerable  volume  of  air  is  left  in  a  highly 
charged  state  as  a  result  of  the  above  processes. 

This  charged  state  will  be  maintained  by  the 
large  potential  gradient  which  energized  tlu 


th rc at c  n i ng  stroke . 
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(4)  The  charged  component  of  air  will  move  in  accordance 
with  the  laws  governing  electrostatic  behavior. 

Positive  ions  will  be  attracted  to  a  negatively 
charged  region  and  vice  versa.  Ion-air  molecule 
collisions  will  set  uncharged  air  in  motion. 

(5)  The  air  set  in  motion  will  have  a  tendency  to  move 
in  a  counterclockwise  vortex  due  to  the  earth's 
magnetic  field  and  the  positive  polarity  of  the  ions. 

(6)  Vortex  motion  will  develop  until  the  aerodynamic 
factors  assume  control. 

(7)  The  tornado  will  end  when  the  local  potential 
energy  region  of  the  storm  has  been  dissipated 

or  the  vortex  has  reached  a  size  of  such  proportions 
that  the  vortical  motion  is  unstable. 

Lewellen  (Ref.  26)  postulates  that  a  magnetohydrodynamic  body 
force  can  be  used  to  drive  a  two-dimensional,  axisymmetric ,  steady 
vortex  in  a  viscous,  conducting  fluid.  In  order  to  drive  such  a 
vortex,  Lewellen  determines  the  field  geometries  which  lead  to 
an  accelerating  body  force  in  the  tangential  direction.  A  force 
of  this  type  results  from  either  an  axial  electric  field  and  radial 
magnetic  field  or  an  axial  magnetic  field  and  radial  electric  field. 

The  theories  discussed  thus  far  postulate  that  the  electrical 
activity  is  a  cause  or  a  driving  force  rather  than  an  effect  of  the 
interactions  of  charged  particles  and  vortices.  Several  authors, 
including  Freier  (Ref.  17)  and  Crozier  (Ref.  18),  Silberg  (Ref.  15), 
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and  Loeb  (Ref.  27)  have  postulated  that  the  electrical  activity 
associated  with  tornadoes,  dust  devils  and  other  vortex  systems  may 
be  the  result  of  vortex  motion  rather  than  the  cause  of  it. 

Freier  and  Crozier,  who  both  measured  the  electric  field 
associated  with  dust  devils,  compared  their  experimental  results  to 
theoretical  results  based  on  a  simple  dipole  theory.  Both  assumed 
a  straight  line  passage  of  a  simple  vertical  dipole  at  a  known 
constant  velocity.  Crozier  obtains  good  agreement  between  the 
experimental  data  and  the  theoretical  curve.  The  simplest  model  of 
a  vertical  dipole  having  the  apparent  electric  moment  of  the  dust 
devil  would  be  a  uniform  distribution  of  negative  space  charge  density 
in  an  approximately  cylindrical  dust  cloud.  Crozier  estimates  that 
the  space  charge  density  of  the  cloud  is  about  1.6  x  10^  elementary 
charges  per  cubic  centimeter.  Although  this  density  may  seem  high 
in  view  of  the  atmospheric  charge  densities  usually  encountered, 

Crozier  postulates  that  the  dust  particle  density  in  such  a  cloud 
may  be  as  high  as  several  thousand  particles  per  cubic  centimeter. 

Since  the  dust  has  just  been  picked  up  from  the  ground  by  a  relatively 
strong  wind,  the  average  particle  size  is  probably  quite  large  and  the 
possibility  exists  that  each  particle  can  carry  several  hundred  elementary 
charges  per  particle.  Thus,  the  estimated  space  charge  densities  in 
the  dust  cloud  arc  not  unreasonable. 

Crozier  suggests  that  a  principal  feature  of  the  charging 
mechanism  is  asymmetrical  electrification  of  the  dust  particles  at 
the  ground  surface,  at  the  initial  separation  from  the  ground  or  on 


17 


subsequent  impacts.  Collisions  between  large  and  small  particles  with 
in  the  cloud  may  also  be  a  factor  in  the  charging  mechanism. 

Crozier  also  attempts  to  explain  the  W- shape  in  the  potential 
gradient  record  of  the  Sahara  Desert  dust  devil  reported  by 
Freier  (Ref.  17).  He  suggests  first  that  it  is  possible  for  the 
electrical  activity  to  increase  rapidly  due  to  rapid  increases  in 
dust  density,  for  example.  Freier  pointed  out  that  the  W-shape  of 
the  curve  could  be  explained  by  assuming  that  there  was  positive 
charge  in  the  lower  portion  of  the  column  and  negative  charge  in 
the  upper. 

In  his  discussion  of  the  dehydration  and  burning  produced  by  a 
tornado,  Silberg  (Ref.  15)  postulates  that  a  pulsed  oscillatory  ring 
current  or  ring  discharge  could  exist  within  a  tornado  cloud  and 
possess  a  radiation  field  that  could  be  responsible  for  these  heating 
effects  at  close  distances  and  the  observed  and  measured  tornado 
sferics  at  far  distances.  The  ring  current  is  defined  as  an  electrode 
less  electrical  discharge  where  the  electric  field  within  the  gas 
forms  closed  curves.  Silberg  cites  several  accounts  of  tornado 
phenomena  which  could  be  very  similar  to  the  induced  ring  discharge. 

He  also  cites  measurements  of  the  radiation  produced  by  the  electrical 
discharges  within  a  tornado  known  as  sferics. 

For  simplicity,  Silberg  assumes  the  ring  current  is  uniform 
and  it  is  sufficient  to  deal  with  a  single  frequency.  He  then 
estimates  the  relatively  low  frequency  radiation  field  at  the  ground 
plane  by  assuming  the  ring  current  behaves  as  a  low  frequency  loop 


antenna.  Silberg  concludes  that  by  using  realistic  values  for  the 
ring  current,  it  is  possible  to  produce  both  the  observed  heating 
effects  at  short  distances  and  the  sferics  measurements  at  larger 
distances  without  unreasonable  electric  field  strengths  or  power 
densities . 

Loeb  (Ref.  27)  postulates  that  the  luminous  phenomena 

occasionally  seen  in  tornadoes  and  the  electrical  charging  of  dust 

devils  can  be  explained  by  considering  the  charging  mechanisms 

existing.  He  cites  the  work  of  Lavan  and  Fejer  (Ref.  22)  who  observed 

that  swirling  jets  of  water  or  saturated  air  enclosed  in  transparent 

insulating  walls  produced  extensive  glow  discharges.  These  discharges 

were  shown  to  be  caused  by  Lenard  spray  electrification.  Larger  drops 

of  water  were  thrown  to  the  insulating  walls  with  the  fine  mist  of 

negatively  charged  droplets  in  the  interior  of  the  tube  as  a  result 

of  centrifugal  forces.  The  potential  gradients  existing  were 

sufficient  to  set  up  a  glow  discharge  in  the  moist  air.  The 

potentials  measured  were  over  20  kilovolts  and  currents  of  90 

microamperes  were  recorded.  Radio  frequency  radiations  of  from 
A  7 

10  to  10  cycles  per  second  were  also  detected  from  the  positive 
corona  streamers.  The  luminosity  arises  from  positive  streamers 


coming  from  spray  disintegrated  drops  in  a  high  electric  field. 
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Loeb  states: 

"That  such  phenomena  may  occasionally  appear  in  the 
tornado  funnel  is  not  surprising  and  depends  on  the 
amount  of  moisture  present  and  entrained.  The  centrifugal 
forces  in  any  case  will  hurl  the  larger  solid  masses,  in¬ 
cluding  positive  water  drops,  toward  the  outside  of  the 
funnel,  and  the  finer  materials,  including  the  negative 
water  spray,  will  be  interior.  If,  in  addition,  solid 
elastic  impacts  of  water  are  needed,  these  would  result 
from  impact  with  larger  debris  particles.  VJhether  a 
particular  funnel  is  luminous  or  not  depends  on  the  water 
entrained--that  is,  on  the  character  of  the  ground 
traversed ." 

He  also  suggests  that  highly  electrically  charged  analogs  of 
tornado  funnels  also  exist  even  without  water  and  cites  the  work  of 
Crozier  (Ref.  18).  Loeb  suggests  that  asymmetrical  contact  charging 
of  dust  with  heavier  particles  is  responsible  for  the  static 
electrification  of  dust  devils. 

The  theories  presented  thus  far  might  lead  one  to  believe  that 
there  is  an  absence  of  hydrodynamic  theories  for  tornadoes  and 
dust  devils.  This  would,  of  course,  be  an  erroneous  conclusion; 
many  hydrodynamic  theories  exist  for  tornadoes,  dust  devils,  and 
other  vortex  systems. 

Turner  and  Lilly  (Ref.  28)  describe  a  carbonated-water  tornado 
vortex  whose  method  of  producing  bouyancy  by  changes  of  phase  in 
the  fluid  itself  is  closely  analogous  to  the  release  of  buoyancy 
by  condensation  in  the  atmosphere.  When  considering  the 
implications  of  their  experimental  results  to  the  atmosphere,  they 
conclude  that  the  most  important  condition  arrears  to  be  the 
presence  of  a  convcctively  overturning  region  above  which  is 
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consistent  with  observations  of  all  tornadoes  and  most  water-spouts 
and  current  theoretical  descriptions  of  tornado  vortices.  Most 
current  theoretical  descriptions  of  tornadoes  agree  that  the 
formation  is  related  to  the  concentration  of  pre-existing  angular 
momentum  by  convective  processes.  The  presence  of  unstable  air 
near  the  ground,  and  condensation  in  the  vortex  probably  contribute 
to  the  persistence  and  intensity  of  the  vortex  but  it  is  possible 
to  form  one  without  it.  It  is  suggested  that  dust  devils  are  the 
result  of  convection  solely  from  an  unstable  region  near  the  ground. 

The  authors  conclude  that  the  mechanism  of  tornado  formation 
appears  to  be  a  concentration  of  the  mean  angular  momentum  by  an 
inflow  to  the  center  of  the  vortex  near  the  bottom  boundary  as  a 
result  of  a  decrease  in  pressure.  This  results  in  an  increased 
velocity  of  circulation  in  the  center  and  at  the  same  time 
probably  a  decreased  circulation  in  a  region  above  the  vigorous 
convection. 

Turner  (Ref.  29)  in  1966  produced  a  laboratory  model  of  a 
tornado  which  incorporates  two  features  believed  to  be  important 
to  the  understanding  of  the  tornado.  First,  he  showed  that  a 
vortex  can  be  driven  from  above  by  a  mechanism  analagous  to 
convection  in  a  cloud,  and  that  density  differences  within  the  funnel 
itself  are  not  essential.  Associated  with  this  mechanism  of 
formation  is  a  vertical  circulation,  with  upflow  in  the  center 
surrounded  by  a  compensating  annular  downflow.  Second,  he  showed 
that  the  bottom  boundary  has  a  strong  influence  on  the  vortex, 
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since  Che  horizontal  and  vertical  flows  are  linked  there  by  a  rapid 
radial  inflow  in  a  thin  boundary  layer. 

'Turner  also  proposed  an  approximate  theoretical  description 
of  the  vortex. 

Logan  (Ref.  30)  developed  a  simple  analytical  model  for 
describing  the  velocity  fields  in  a  dust  devil.  The  model  has 
a  viscous  inner  region  or  boundary  layer  composed  of  a  Prandtl 
layer,  a  layer  in  which  there  is  a  transition  from  a  viscous  boundary 
layer  to  an  inviscid  outer  region  and  an  inviscid  "outer"  region 
which  observations  indicate  can  be  approximated  by  a  Rankine  com¬ 
bined  vortex.  Logan's  model  cannot  properly  describe  the  core 
downdraft  and  does  not  deal  explicitly  with  the  thin,  unstable, 
thermal  boundary  layer.  He  suggests  that  a  more  sophisticated  model, 
including  an  energy  equation  and  equation  of  state,  would  yield 
temperature  and  density  fields. 

Barcilon  (Ref.  31)  proposes  a  theoretical  and  experimental  model 
for  a  dust  devil.  He  considers  the  flow  field  found  in  a  steady, 
axially  symmetric,  strong  atmospheric  vortex  of  the  dust  devil  type. 
The  plane  containing  the  axis  of  the  vortex  is  divided  into  an 
inviscid  region  where  an  unstratified  free  vortex  flow  exists,  a 
plate  boundary  layer  region,  a  corner  region,  and  an  axial 
boundary  layer  region.  Each  region  is  considered  in  turn  and  a 
solution  to  the  overall  flow  field  is  found  by  matching  the  flow 
field  at  the  interface  between  two  adjacent  regions.  Using  this 
model,  Barcilon  predicted  some  of  the  features  found  in  the  model 
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where  the  inviscid  flow  is  stratified.  He  then  discusses  an  experi¬ 
mental  demonstration  in  which  the  fluid  is  given  an  unstable 
temperature  stratification  with  height  and  a  source  of  angular 
momentum.  Under  certain  conditions  a  vortex  of  the  dust  devil  type 
formed . 

Battan  (Ref.  32)  notes  that  a  dust  devil,  unlike  a  tornado, 
forms  in  the  absence  of  clouds  and  as  a  consequence  derives  none 
of  its  energy  from  latent  heat.  Calculations  have  been  made  by 
him  to  determine  the  total  amount  cf  energy  (internal  and  potential) 
available  for  conversion  into  kinetic  energy,  and  the  amount  of 
kinetic  energy  in  a  dust  devil.  He  concludes  that  his  calculations 
reaffirm  the  fact  that  the  kinetic  energy  of  a  dust  devil  results 
from  a  reduction  of  the  internal  and  potential  energy  of  the  air 
involved . 

Dergarabedian  and  Fendell  (Ref.  33)  discuss  the  parameters 
governing  the  generation  of  free  vortices.  They  propose  an 
expansion  of  the  Navier-Stokes  equations  which  provides  a  new 
basis  for  free  vortex  equations.  The  description  of  vortex 
generation  leads  to  two  main  parameters:  the  first  is  the  Eknan 
number,  E  =  v/F^,  where  v  is  the  kinematic  viscosity  and  is  the 
ambient  circulation;  the  second  indicates  that  unless  the  product 
of  the  area  of  the  updraft  and  the  average  vertical  gradient  of 
the  vertical  velocity  greatly  exceeds  the  kinematic  viscosity 
coefficient,  no  localized  vortical  flow  develops.  The  model  is 
for  the  most  part  axially  invariant,  but  the  lines  of  constant 


pressure  described  by  the  model  are  shown  to  have  a  funnel-cloud¬ 
like  shape.  The  authors  propose  that  the  analysis  may  describe 
the  competition  between  the  convective  influx  to  concentrate  the 
ambient  circulation  about  the  axis  of  the  updraft  and  the  efforts 
of  the  viscous  forces  to  smooth  out  the  circulation  distribution. 

Donaldson  and  Sullivan  (Ref.  34)  made  a  study  of  the  complete 
class  of  solutions  of  the  Navier-Stokes  equations  in  which  the 
radial,  tangential  and  axial  velocities  in  cylindrical  coordinates 
are  of  the  form  u  =  u(r),  v  =  v(r),  and  w  =  zw(r).  These  solutions 
were  found  to  represent  a  large  class  of  three  dimensional  viscous 
vortex  motions.  The  solutions  obtained  show  that  vortex  motions 
are  possible  which  have  more  than  one  cell.  That  is,  the  flow  may 
not  simply  spiral  in  toward  an  axis  and  out  along  it  as  in  a  one- 
celled  configuration  but  may  have  nested  regions  of  successively 
reversed  axial  flow.  The  flow  streamlines  and  velocity  profiles  for 
vortices  having  one  and  two  cells  are  shown  in  Figure  1.  The 
behavior  of  the  solutions  in  passing  from  single  to  multiple  celled 
configurations  is  discussed  and  the  solution  for  the  case  of  a  two- 
celled  analogue  to  Burgers'  unconstrained  vortex  (Ref.  35)  reported 
by  Sullivan  (Ref.  36)  is  given.  A  solution  of  the  type  reported  by 
Sullivan  is  of  considerable  practical  importance  because  many 
laboratory  vortex  chambers  exhibit  a  two-celled  flow  and  many  large 
vortex  phenomena  such  as  hurricanes  and  tornadoes  also  exhibit  this 
two-celled  character. 

The  discussion  now  turns  to  look  at  some  of  the  related 
electrical  phenomena,  such  as  discharge  stabilization  and  charging 
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mechanisms,  which  should  give  better  insight  into  the  problem  of 
determining  the  role  of  electrical  activity  in  vortex  systems. 

Vonnegut,  Moore,  and  Harris  (Ref.  37)  conducted  simple  laboratory 
experiments  to  determine  what  effects,  if  any,  a  vortex  would  have 
on  a  continuous  high  voltage  discharge.  They  found  that  the  vortex 
produced  the  following  effects  on  the  high-voltage  discharge: 

1)  it  makes  the  discharge  more  stable; 

2)  it  causes  the  character  of  the  discharge  to  change 
from  successive  sparks  to  a  continuous  glow; 

3)  it  increases  the  distance  over  which  the  discharge 
can  be  maintained; 

4)  it  decreases  the  electrical  resistance  of  the  discharge; 

5)  it  decreases  the  noise  of  the  discharge; 

6)  it  decreases  the  production  of  the  electromagnetic 
radiation  in  the  discharge  band; 

7)  it  causes  the  discharge  to  center  itself  along  the 
axis  of  rotation. 

The  authors  believe  that  the  effects  produced  by  the  vortex 
on  the  discharge  can  be  explained  in  terms  of  the  centrifugal  force 
it  exerts.  A  radial  pressure  gradient  which  is  zero  at  the  axis 
and  increases  with  the  distance  away  from  it  is  produced  by  the 
rotation  of  the  air  in  the  tube.  In  a  pressure  gradient  like  this, 
bouyant  forces  cause  the  hot  gas  in  the  discharge  to  move  toward 
the  region  of  lowest  pressure  at  the  axis.  The  centrifugal  action 
automatically  confines  the  discharge  at  the  center  of  the  vortex. 

This  stabilizes  the  vortex  by  forcing  the  hot  ionized  gases  together 
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at  the  axis  and  by  preventing  them  from  mixing  with  the  cool 
surrounding  air. 

The  authors  also  concluded  that  the  results  suggest  that  the 
tornado  funnel  might  serve  as  a  preferred  path  for  lightning  and  that 
the  vortex  may  stabilize  and  modify  lightning  discharges. 

Wilkins  and  McConnell  (Ref.  38)  in  1968  attempted  to  determine 
experimentally  the  threshold  conditions  for  vortex-stabilized 
discharges  in  the  atmosphere  and  whether  or  not  a  vortex  core  is 
an  improved  conductor  for  arc  discharges.  They  concluded  that  if 
point  discharging  is  a  factor  in  the  prevention  of  lightning  as 
indicated  by  the  laboratory  analogy,  then,  since  high  velocity 
flow  tends  to  speed  up  point  discharging,  a  tornado  funnel  is  not 
a  favored  region  in  this  respect  for  cloud  to  ground  lightning. 

The  core  region  is  however,  more  favorable  than  the  region  just 
outside  it,  and  if  the  air  above  the  tornado  were  stressed  to  the 
breakdown  point,  it  is  possible  for  a  ground  stroke  to  occur  along 
the  vortex  axis. 

The  experiments  show  that  the  threshold  power  required  to 
maintain  a  steady  discharge  decreases  as  vortex  velocity  increases, 
but  this  effect  is  not  very  pronounced  and  one  would  not  conclude 
from  it  that  a  tornado  axis  would  provide  a  highly  favored  environment 
for  continuous  cloud- to-ground  discharging.  The  power  for  the  initial 
lightning  discharge  would  somehow  have  to  be  maintained  for  periods 
of  several  seconds,  and  there  is  no  reason  to  suppose  this  is  possible. 
The  experiments  indicate  that  the  discharge  will  flame  out  between 
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lightning  strokes  because  the  centrifuge  action  of  the  vortex 
in  confining  and  concentrating  the  plasma  is  not  sufficient  to 
appreciably  lengthen  the  restrike  time.  In  view  of  these  considerations 
the  authors  suggest  that  a  tornado  funnel  is  not  a  highly  favored 
region  for  multistroke  lightning  discharges  unless  the  tornado 
thunderstorm  is  a  considerably  more  efficient  charging  system 
than  other  thunderstorms. 

Dunham  (Ref.  39)  made  measurements  of  the  charge  produced  by 
snow  striking  various  surfaces  at  high  speed.  He  concluded  that 
electrostatic  charging  by  friction  and  fragmentation  of  snow 
particles  is  dependent  on  the  substance  the  snow  is  striking  and 
the  charging  rate  increases  with  a  decrease  in  temperature  from 
0°C  to  -7°C.  Charging  rate  per  gram  varies  as  the  square  of  the 
speed. 

Latham  (Ref.  AO)  made  an  analysis  of  the  available  observations 
on  the  electrical  properties  of  natural  and  artificial  snowstorms. 

These  electrical  effects  are  in  qualitative  agreement  with  the  theory 
of  charge  transfer  associated  with  temperature  gradients  in  ice, 

*  i 

which  has  been  shown  to  provide  a  tenable  theory  of  thunderstorm 
electrification.  The  electrical  effects  associated  with  natural  and 
artificial  sandstorms  have  a  marked  resemblance  to  the  snowstorm 
electrification  phenomena  and  appear  to  be  explainable  in  terms  of  a 
similar  temperature  gradient  effect. 

Kunkel(Ref.  Al)  conducted  an  extensive  study  on  the  static 
electrification  of  dust  particles  on  dispersion  into  a  cloud  under 
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a  variety  of  conditions  ranging  from  blowing  with  a  minimum  of 
turbulence  to  violent  and  maximum  impact  with  various  types  of 
surfaces.  He  concluded  that  there  appears  to  be  no  doubt  that  the 
act  of  separating  a  mass  of  solid  particles  in  contact  with  each 
other  into  a  cloud  of  particles  of  assorted  sizes  yields  in  all  cases 
a  cloud  which  as  a  whole  is  neutral,  but  in  which  practically  all  the 
particles  are  charged.  The  number  of  particles  of  similar  size  but 
with  opposite  charges  was  about  equal.  The  magnitude  of  the  charges 
did  not  increase  as  rapidly  as  the  surface  of  the  particle.  No 
correlation  was  observed  between  particle  size  and  sign  of  charge. 
There  appears  to  be  strong  evidence  that  charging  occurs  on  separation 
of  the  contacts  between  particles  in  the  dispersion  of  the  cloud. 
Charging  did  not  appear  to  be  affected  by  humidity.  Heterogeneous 
systems  making  contact  with  solid  walls  of  different  composition 
from  the  powder  gave  consistent  asymmetry  of  charge  of  varying  degrees 
depending  on  the  proportion  of  particles  striking  the  surfaces 
relative  to  those  just  separated. 

Much  work  has  been  done  in  the  area  of  raindrop  electrification 
and  raindrop  charge  and  thunderstorm  electrification  by 
Gunn  (Ref.  42,  43,  44,  45,  46).  His  work  is  significant  to  this  study 
in  that  the  droplet  charging  mechanisms  he  presents  and  discusses  are 
possibly  the  mechanisms  causing  the  initial  electrification  of 
droplets  in  a  thunderstorm  and  buildup  of  the  potential  gradients 
in  thunderstorms. 

It  is  hoped  that  the  accounts  of  electrical  activity  in  vortex 
systems,  the  electrical  and  hydrodynamic  theories  for  tornadoes, 
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dust  devils,  and  vortex  systems,  and  the  charging  mechanisms  related 
to  other  electrical  phenomena  presented  here  will  permit  some 
substantial  conclusions  to  be  drawn  concerning  the  role  of  electrical 
activity  in  vortex  systems. 

The  first  area  to  be  discussed  is  the  role  of  electrical  activity 
in  tornadoes.  It  appears  that  three  choices  are  available  here: 

1)  tornadoes  are  caused  by  electrical  activity;  2)  electrical 
activity  is  caused  by  tornadoes;  3)  any  association  between  the  two 
is  fortuitous.  Each  of  these  choices  will  be  examined  in  greater 
detail.  If  tornadoes  are  caused  by  electrical  activity,  than  the 
electrical  heating  theories  of  Vonnegut  (Ref.  9)  and  Colgate  (Ref.  24), 
the  electromagnetic  theory  of  Rathbun  (Ref.  25),  the  idea  of  electro¬ 
dynamic  acceleration  of  ions  proposed  by  Vonnegut  (Ref.  9),  or  the 
magnetohydrodynamically  driven  vortex  theory  of  Lewellen  (Ref.  26) 
could  be  used  to  explain  the  relationship.  The  electrical  heating 
theory  does  not  appear  feasible  because  it  requires  long,  sustained 
lightning  strokes  or  some  totally  unknown  phenomenon  to  heat  the  air 
sufficiently  to  provide  the  intense  winds  fouhd  in  a  tornado.  A 
tornado  funnel  may  be  a  preferred  path  for  a  lightning  stroke,  but 
it  is  doubtful  that  even  this  could  lead  to  a  very  large  heating  of 
air.  This  theory  is  certainly  not  sufficient  to  account  for  the 
intense  updrafts  and  high  angular  momentum  associated  with  funnel 
clouds . 

Rathbun' s  (Rof.  25)  electromagnetic  theory  suffers  because  there 
is  no  known  physical  process  that  could  cause  the  ion  concentration 
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in  the  air  to  remain  as  large  as  the  theory  assumes  it  is  for  as 
long  a  period  of  time  as  the  theory  assumes  also.  The  idea  of 
electrodynamic  acceleration  of  ions  or  a  magnetohydrodynamically 
driven  vortex  are  both  dependent  on  the  development  of  a  tangential 
body  force  within  the  air  to  accelerate  it.  Lewellen  (Ref.  26)  has 
shown  that  this  requires  an  axial  electric  field  and  a  radial 
magnetic  field  or  conversely,  a  radial  electric  field  and  an  axial 
magnetic  field.  The  first  combination  can  be  eliminated  because  as 
Wilkins  (Ref.  23)  has  pointed  out,  there  is  no  way  for  a  radial 
magnetic  field  to  exist.  The  second  combination  is  physically 
possible,  but  it  is  felt  that  the  axial  electric  field  required  to 
maintain  the  ion  concentration  will  probably  overshadow  any  radial 
elec  trie  field. 

With  the  theories  available  at  present,  it  appears  unlikely 
that  tornadoes  are  definitely  and  directly  caused  by  electrical 
activity.  It  may  be  possible  that  electrical  forces  help  to  maintain 
a  tornado,  but  it  is  difficult  to  conclude  this  based  on  present 
day  knowledge  or  thinking. 

The  third  hypothesis:  any  association  between  tornadoes  and 
electrical  activity  is  fortuitous,  will  be  discussed  next.  The 
various  accounts  of  electrical  activity  in  tornadoes  make  it 
difficult  to  conclude  that  this  activity  happens  by  chance.  Present 
knowledge  in  the  field  of  thunderstorm  electrification,  raindrop 
electrification  and  other  atmospheric  electrical  phenomena  as  evi¬ 
denced  by  the  work  of  Gunn  (Ref.  A  6) ,  for  example,  quickly  dispel 
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any  notions  that  the  electrical  activity  is  fortuitous.  The  relation¬ 
ship  may  be  a  complex  one  with  many  interacting  variables,  but  there 
is  no  doubt  that  a  relationship  does  exist. 

The  second  hypothesis:  the  tornado  or  vortex  motion  causes 
electrical  activity,  appears  to  be  most  plausible  at  this  time.  The 
remarks  by  Loeb  (Ref.  27),  the  work  of  Lavan  and  Fejer  (Ref.  22),  the 
measurements  by  Freier  (Ref.  17),  Crozier  (Ref.  18),  and  Bradley  and 
Semonin  (Ref.  19)  all  serve  to  bolster  the  belief  that  this  is  the 
correct  relationship  between  the  two  phenomena.  The  measurements  of 
Freier  (Ref.  17),  Crozier  (Ref.  18),  and  Bradley  and  Semonin  (Ref.  19) 
on  the  electric  fields  associated  with  dust  devils  all  indicated  that 
this  vortex  system  had  a  large  region  of  negative  charge  above  it. 

Loeb  (Ref.  27)  credits  this  to  asymmetrical  contact  charging  of  dust 
with  heavier  particles.  The  centrifugal  action  of  the  vortex  and 
gravitational  attraction  would  cause  the  lighter,  negatively  charged 
dust  to  be  at  the  top  of  the  dust  devil.  The  work  of  Lavan  and 
Fejer  (Ref.  22)  on  swirling  flows  in  ducts  showed  that  swirling  jets 
of  saturated  air  or  water  enclosed  in  insulating  walls  produced  glow 
discharges.  Here  the  centrifugal  action  of  the  vortex  served  as  a 
charge  separator  by  causing  the  larger,  positive  water  drops  to  go  to  the 
wall  with  a  fine  mist  of  light  negative  charges  remaining  in  the  center 
of  the  tube.  A  glow  discharge  in  the  tube  was  set  up  by  the  more  than 
adequate  potential  gradients  that  developed. 

This  survey  has  served  to  illustrate  that  the  determination  of  the 
role  of  electrical  activity  in  vortex  systems  such  as  the  tornado  or 
the  dust  devil  is  a  very  complex  problem,  requiring  esearch  in  many 
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different  areas.  This  study  of  charged  particle  motion  in  a  free- 
vortex  flow  field  is  intended  to  determine  if  particle  motion  in  a 
free  vortex  can  be  influenced  by  the  interaction  between  the  particle 
charge  and  an  applied  electrostatic  field. 


DISCUSSION  AND  DESCRIPTION  OF  PARTICLE-FLOW  SYSTEMS  AND  THEORY 


In  order  to  simplify  the  analysis  of  the  motion  of  a  charged 
particle  in  a  free-vortex  under  the  influence  of  an  applied 
electrostatic  field,  several  assumptions  were  made.  First,  rather 
than  attempting  to  analyze  the  motion  of  a  large  number  of  particles 
having  some  size  distribution  and  charge  distribution,  this  analysis 
is  concerned  with  the  motion  of  a  single  charged  particle.  The 
analysis  of  the  general  motion  of  a  gas-solid  system  with  some  size 
distribution  of  solid  particles  and  including  interactions  among 
particles  may  be  treated  using  the  methods  discussed  by  Soo  (Ref.  47). 
Second,  the  particle  is  assumed  to  be  charged  prior  to  its  intro¬ 
duction  into  the  free-vortex  flow  field.  No  explanation  of  how 
the  particle  became  charged  will  be  advanced,  other  than  to  indicate 
that  electrification  of  solid  particles  does  occur  when  contact  and 
separation  are  made  between  the  particles  and  a  dissimilar 
material,  or  a  similar  material  under  different  surface  conditions. 
Kunkel  (Ref.  41)  has  made  an  extensive  study  on  the  static  electri¬ 
fication  of  dust  particles  on  dispersion  in  a  cloud  in  which  he 
concludes  that  electrification  occurs  when  contact  and  separation 
are  made  between  the  particles.  Third,  the  assumption  is  made  that 
the  flow  field  is  an  axisymmetric  free-vortex,  that  is,  the  tangential 
flow  is  inversely  proportional  to  the  radius.  To  avoid  the  problem 
of  the  infinite  velocity  at  the  origin,  the  particle  is  assumed  to  be 
initially  at  some  finite  radius  and  its  predominant  motion  is  away 


from  the  origin.  Fourth,  the  assumption  is  made  that  the  drag  force 
on  the  particle  is  proportional  to  its  relative  velocity  and  that 
the  drag  coefficient  can  be  represented  using  Stokes'  law.  This 
assumption  is  valid  at  low  Reynolds  number  when  streamline  or  viscous 
flow  exists  around  the  particle.  The  drag  force  is  given  by: 


CDP 


tt  -  v2 
2 
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and  the  drag  coefficient  is  given  by: 


=  24u 

D  "  p  |  u  -  v  |  D 
P1  P 


Therefore,  the  drag  force  can  be  expressed  as: 


Fd  =  3WDpfl  -  vp)  . 

Fifth,  it  was  assumed  that  the  applied  electrostatic  field  was  the 
result  of  a  line  source  of  positively  charged  particles  located  at 
the  axis  of  the  free-vortex.  The  electrostatic  field  due  to  the 
charged  line  source  is  given  by  E  =  X/2nc^r  (Ref.  48). 

In  order  to  analyze  the  motion  of  a  charged  particle  in  a  free- 
vortex  under  the  influence  of  an  applied  electric  field  and  determine 
the  parameters  affecting  the  particle  motion,  four  different  particle 
flow  systems  were  considered.  The  first  case  analyzed  the  motion  of 
uncharged  particle  in  a  free-vortex.  This  system  was  analyzed  to 
provide  a  baseline  for  determining  the  effects  of  the  drag  force  and 
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centrifugal  force  on  particle  motion.  The  analysis  of  this  system 
is  referred  to  as  Case  1  in  the  following  discussions.  The  second 
system  considered  consisted  of  a  positively  charged  particle  free  to 
move  in  a  viscous  medium  under  the  influence  of  an  applied  electro¬ 
static  field  due  to  a  line  source  of  positively  charged  particles. 

The  analysis  of  this  system  indicated  the  effects  of  the  electrostatic 
body  force  on  particle  motion.  The  analysis  of  this  system  will  be 
referred  to  as  Case  II.  The  third  system  investigated  consisted  of 
a  positively  charged  particle  moving  in  a  free-vortex  under  the 
influence  of  an  applied  electrostatic  field.  In  effect,  this  system 
is  created  by  superimposing  the  first  two  systems  considered  on  each 
other.  The  analysis  of  this  system  is  referred  to  as  Case  III,  A 
variation  to  this  system,  in  which  the  particle  is  assumed  to  be 
negatively  charged,  will  be  referred  to  as  Case  III  A.  Finally  the 
fourth  particle-flow  system  consisted  of  two  charged  particles  moving 
in  a  free-vortex  under  the  influence  of  the  electrostatic  field  due 
to  the  charge  on  each  of  the  particles.  The  analysis  of  this  system  is 
referred  to  as  Case  IV. 

The  determination  of  the  equations  governing  the  particle 
motion  for  each  of  the  four  systems  considered  will  now  be  presented. 
The  nomenclature  used  here  is  defined  in  Appendix  A. 
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Case  1 

The  prediction  of  the  motion  of  uncharged  particles  in  a  free- 
vortex  system  follows  the  analysis  of  Hirschkron  and  Ehrich,  (Ref.  1), 
who  studied  the  trajectories  of  uncharged  particles  entrained  in  a 
swirling  flow.  The  following  assumptions  were  made  in  their  analysis: 

1)  The  flow  field  is  an  axisymmetric  free-vortex  with  a  constant 
axial  velocity; 

2)  The  particle  mass  flow  is  small  compared  to  the  carrier- 
fluid  mass  flow; 

3)  The  particles  do  not  interact  with  each  other,  nor  do  they 
affect  the  flow  field; 

4)  The  particles  are  not  charged; 

5)  The  drag  force  is  given  by  the  Stokes  relation: 

-  3  TT/X.  "Dp  (  U.  -  Vp  )  . 

The  equation  of  motion  of  a  particle  is  then  given  by: 

f*  ^  =  3tt/aDp  ( U.  -  '/p')  +  B,  (1) 

where  B  is  the  body  force. 

The  particle  coordinates,  particle  velocity  components,  and  the 
flow  field  components  in  the  (r,Q)  plane  and  the  (r,z)  plane  are 
shown  in  Figures  2  and  3. 

In  (r,6,z)  coordinates,  Equation  (1)  can  be  expanded  to, 

s  (2) 

<»'r(Vr  r  3w.«ty  (•«,-  V,)  }  (3) 
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-  3tt/aDp  (  'O  ~  rtv* p  J  _ 

(4) 

Assume 

that  the  flow  field  is  an  axisymmetric  free-vortex. 

Therefore, 

Mr*  o, 

(5a) 

^  Q  ~  Cv  / 

(5b) 

((.^  =  Ugei  constant. 

(5c) 

Define 

C  5  3tT/4  J^p  ( C(t\  p  . 

(6) 

Applying  5a,  5b,  5c,  and  6  to  Equations  2,  3,  and  A  reduce  the 


latter  to  the  following  form: 

=  -CV,,  O) 

Vr  +  Vc^S.  =  £  ( Cv/r  -  \4Y  (8) 

dr  r  1 

V,  =  C(U?.-V,)  -  J.  (9) 


Setting 

Xq  -  “  ^v/r, 

(10) 

and  substituting  for  vq  in  Equation  (8),  the  following  result  is 

obtaine  d : 

• 

Vr  4  Vr/„  =  -  (?  Xe  . 

Jr  r  e  5 

(11) 

or, 

Vr  J  (hX^')  =  -  C  Xq  . 
r  Jr 

(12) 
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The  variable  (rxg)  can  now  be  separated  such  that 


(k(v  Xq)  _  _  d\T  t  (1“*) 

( rx,")  ~  Vr 

Along  the  particle  trajectory,  VV  =  ^  *  (14) 

Therefore,  ^  l =  —  C  .  (15) 

(rXe'j 

Integrating  Equation  (15)  from  t  -  0  to  t  ■  t, 

,  V  -Ct 

(rx0)  --  ,  ae) 

where  (f  Xe)0  s  (  V«#  -  (17) 

Since  Vs  *  Xa  +  £v/V, 


/e  *  [r,fV.s-  &/i:)e’,t+  <v]/ r.  (18) 

Substituting  for  vQ  in  Equation  (7),  and  applying  Equation  (14), 

•  t  *1^ 

d  \ r  _[  >o(V«.  -f  Cy  |  _  _  /)  v  _  (19) 

r‘  "  ~  u  dt  ' 


Along  a  particle  trajectory,  \/a  - 

0  dt 

Therefore,  =  fo(  ^  '  Cr/lOg"  *.*■<*  . 

dt  r*- 


(20) 

(21) 


Define  v  = 
z 


<(fc2 


dz/dt  and  substitute  for  v 

z 

--  cau.  -  &)  -  ?• 


in  Equation  (4). 


(22) 


Equations  19,  21,  and  22  describe  the  particle  motion  within  the 
free-vortex  system.  The  case  where  the  initial  particle  injection 

velocity,  Vg  ,  is  zero  will  be  examined  in  some  detail, 
o 
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I 


Using  Vq  =  0,  the  equations  describing  the  particle  motion  are 
o  „ 

&JC.  4-  C  I  -  g~  Cj).  ]  =  q 

dt1  dt  r*  > 


(23) 


_  Cv  (i-  e'Ci) .  _ 


dt 


O. 


(24) 


ft  ♦  < &  -  “«•  4r° 


(25) 


These  equations  can  be  nondimensionalized  through  the  use  of 
nondimensional  parameters  R  =  r/rQ,  Z  *  z/zq,  an<^  T  *  Ct. 


+  4E  -  /£vf  (l -  s'7)’’  =  o, 

aT 


ie  _  ^evJO-e  )  = 


dt 


d'z 


o 


+  -  iiifl  +  JL  =  o  , 

j(Tx  dr  Cz.  c*e. 


(26) 


(27) 


(28) 


Assume  the  following  boundary  conditions: 

fit  t  -  0> 

t-B, 

9*  0. 

2  *  ^  "  <f£  ' 

In  nondimensional  form,  the  boundary  conditions  are: 

fit  T-o, 

v-  \  ^  ^  o . 

r%  rtT  > 

^  O, 

.  1  ii  0. 

Equations  26  and  27  cannot  be  explicitly  solved  for  the  particle 
trajectories;  however,  they  can  be  integrated  relatively  easily  using 
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a  numerical  procedure.  A  fourth-order  Runge-Kutta  scheme  was  pro¬ 
grammed  for  a  digital  computer  to  obtain  the  solutions  to  these 
equations.  The  computerized  equations  and  the  Runge-Kutta  scheme 
will  be  described  in  greater  detail  in  the  chapter  on  the  computer 
programs. 

Equation  (28)  can  be  integrated  directly,  subject  to  its 
associated  boundary  conditions. 

£3.  +  4i  -  it*.  +  1:0  (29) 

dr  Cz0  Cfo 

Boundary  Conditions: 

fit  T-  o,  Z  *  1 ,  %  •  O. 

* 

The  solution  to  this  differential  equation,  subject  to  its 
boundary  conditions,  is  given  by: 

2=1+  (^^)(  T  +  e  T- 1)  .  oo) 

In  dimensional  form,  the  axial  particle  trajectory  is  given  by 

H:  +  (C*u-Q)(t/c  +  eCt/c-  l/c)  . 

c 


(31) 
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Case  II 

The  second  case  analyzed  considered  the  notion  of  a  charged 
particle  due  to  the  presence  of  a  line  source  some  distance  rQ 
away  from  the  charged  particle.  Viscous  drag  in  the  gas  is  to  be 
accounted  for  through  the  use  of  the  Stokes  relation.  This  analysis 
is  similar  to  that  conducted  by  Soo  (Ref.  6).  Figure  4  indicates  that 
the  line  source  of  strength  X  coincides  with  the  Z-axis  and  the 
particle  having  charge* to-mass  ratio  (q/m)  is  initially  a  distance 
ro  away  from  the  line  source. 

The  following  assumptions  were  made  in  this  analysis: 

1)  No  flow  field  is  acting; 

2)  The  particle  is  charged; 

3)  The  applied  electrostatic  field  due  to  a  line  source  affects 
the  particle  motion; 

4)  The  drag  force  is  given  by  the  Stokes  relation; 

F*  =  3TT/U  Dp  (  (4  -  tfpV 

The  equation  of  motion  of  the  particle  is  then  given  by: 

mr\p  tdl  -  3  ir/i  £>P  ( tJL  -  7, )  (32) 

Since  the  applied  electrostatic  field  due  to  the  charged  line 
source  acts  in  the  radial  direction  only,  and  since  the  initial 
particle  velocities  in  the  tangential  and  axial  directions  are  zero, 
the  equation  of  motion  reduces  to: 

01  ‘  p  ^  ™  “  “  3n/<^  '4*  +  E''  <33) 


Charged  Line  Source 


,  Charged  Particle 

/  * 
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For  a  line  source,  Er  is  given  by  (Ref.  48): 


Er  = 

(34) 

Defining 

ii 

> 

dr 
dt  5 

(6) 

C  = 

377/1 

(14) 

T?« 

^ka/*re.mr. 

(35) 

and  substituting  Equations  (6),  (14),  (34),  and  (35)  into  Equation  (33), 
tie  following  equation  results: 


dzr  o  dx  D 
L  &+.  “  r 


(36) 


Assume  the  following  boundary  conditions: 

/)t  t  --  o,  rQ> 

This  differential  equation  and  its  associated  boundary  conditions 
can  be  nondimensionalized  through  the  use  of  the  nondimensional 
parameters  R  =  r/r  ,  and  T  =  Ct. 


£p  +  =  0 

c(  T*  <*T  C*R 


(37) 


The  boundary  conditions  are: 

fit  T-o,  K-l,  £*o. 

Equation  (37)  can  be  integrated  using  the  numerical  procedure 


mentioned  in  Case  I. 
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Case  III 

The  third  case  studied  involved  a  charged  particle  moving  in  a 

l  J 

free-vortex  flow  field,  under  the  influence  of  an  applied  electro¬ 
static  field  due  to  an  axial  line  source  of  charges.  The  forces  on 
the  particle  under  consideration  are  the  drag  force,  a  body  force  due 
to  the  applied  electrostatic  field,  and  the  gravitational  force. 

The  following  assumptions  were  made  in  this  case: 

1)  The  flow  field  is  an  axisymmetric  free-vortex  with  a 
constant  axial  velocity; 

2)  The  particle  mass  flow  is  small  compared  to  the  carrier 
fluid  mass  flow; 

3)  The  particle  is  charged; 

4)  The  applied  electrostatic  field  due  to  the  line  source 
affects  the  particle  motion; 

5)  The  drag  force  is  given  by  the  Stokes  relation, 

Fd  -  »Dp  (S  -  Vp). 

Figures  5  and  6  indicate  the  particle  coordinates,  particle  velocity 
components,  and  fluid  velocity  components  in  the  (r,6)  and  (r,z) 
planes.  Also  indicated  in  Figures  5  and  6  are  the  initial  particle 
positions  and  the  location  of  the  charged  line  source  on  the  Z-axis. 
The  equation  of  motion  of  the  particle  is  then  given  by: 

+s  (38) 

£ It 

where  B  is  the  body  force. 

Equation  (38)  can  be  expanded  in  (r,6,z)  coordinates  to 

fi 

JF  -  ^^(Ur-Vr)  +  Ev  O’) 

[] 


FIGURE  5 


Charged  Line  Source 


FIGURE  6 
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*  XX)  =  B3, 

m'f('/'Tr)  =  Sn^-Cpf;/,- O  ,  B, 


(40) 

(41) 


In  the  radial  direction,  Br  =  xqpgc/2neor;  in  the  tangential 
direction,  B0  =  0;  and  in  the  axial  direction,  E  =-m  p 

2  pb " 

Assuming  that  the  flow  field  is  an  axisymmetric  free-vortex, 

ilr  -  O , 

tV  c-/r 


(5a) 

(5b) 


^ i  *  ^2C 


constant 


(5c) 


Substituting  into  Equations  39,  40,  and  41  for  B  i>  b 

r  »  n0  >  z  > 

and  ur.  UQ,  Uz,  the  following  equations  result: 

4  #  -  £  *  C  (-  v.)  +  I-  , 


*  =  0  (  k  -  v,)  , 

where  C  =  3TVD  /m  and  D  =  ). g  q  /2nr  m  y  -  _  . 

P  P  bcV  oV  Let  X0  ve"Gv/r,  and 

substitute  for  Vq  in  Equation  43. 


(42) 

(43) 

(44) 


Vr 


^  IT 


lU, 

r 


-  ex. 


(45: 


Equation  (45)  can  be  rewritten  as: 


X.  J 


<•  ) 


-  v 

/  V  (. 


(46) 


48 


The  variable  rXg  can  now  be  separated  such  that: 

j(y  _  q  d.  v 

(rxe)~  Vr 


Using  =  dr/dt,  d(lrXe')  ^  ^ 

(rx«) 

Integrating  Equation  (48)  from  t  =  0  to  t , 

-  ct 

IT  XQ  =  £  Xe0  > 

where 

r.  r.(  Vs,  -  C,/0. 

Since 

(47) 

(48) 


(49) 


V©  -  v:(Ve0-  a/Oe"ct  + 

r  ■  r 


(50) 


Substituting  for  vQ  in  Equation  (42),  using  =  dr/dt,  the  equation 
for  the  radial  trajectory  becomes: 

c/zr  _  [r.  ( N/e.  -  G/ VQ  4  Cv1\  /,  dr  _ 


O 


r’  '  "  <<t  r 

Along  the  particle  trajectory  Vq  =  r  d9/dt.  The  equation  for  the 
tangential  particle  trajectory  becomes: 


(51) 


"  Cv/c. ) 0  +Cv 

r  It  "  r 

le.  -  ’■£  C  V^a  -  Cv/r)e  ct-r  Cv 


or, 


dt 


I' >- 


=  o. 


(52) 

(53) 


Define  v  =  dz/dt  and  substitute  for  v in  Equation  (44) 


C  ii  -  C  +■  A  =  o . 

dt  c 


.(! 

<al 


(54) 
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Equations  51,  53,  and  54  describe  the  motion  of  a  charged 

particle  in  a  free-vortex  system  under  the  influence  of  an  applied 

electrostatic  field.  Once  again  the  initial  particle  velocity  in  the 

tangential  direction,  Vq  .will  be  assumed  zero.  This  case  will  be 

o 

examined  in  some  detail. 

Using  v0  =  0,  the  equations  of  motion  are: 


i!l  nr  C  ^  - 
d  t*  it 


)-  e 


-  «  o 

r 


de  _  'Cv(\- e’ci\  s  o, 
it 

•  +C«-  CU:  +  ]*0. 


dzi  .  n  6  2 


The  nondimens ional  form  of  these  equations,  using  R  =  r/rQ,  T  =  Ct, 


z/zo, 

is  given  by: 

dTx 

dT 

dQ 

./Cv  \ 

i-e"T 

dT 

\crS  J 

1  K 

JM 

■t  ^ 

(J.Zq 

iiz 

<dT 

C20 

(Croy-  K 


Assume  the  following  boundary  conditions: 

At  t®  C, 

r-  Vr 1  dr/d t  ®-  o, 

O'-  o, 

.  ^  x/j,  -  i~/dt  * 

In  nondimens ional  form,  the  boundary  conditions  are: 
At  T= 

fM>  df./dT - O, 

e-'O, 

2  -  1 }  d  0  • 
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liquations  (58)  and  (59),  describing  the  radial  and  tangential 
particle  trajectories,  can  be  integrated  using  the  numerical 
procedure  mentioned  previously  in  Case  1.  The  computerized  form 
of  these  equations  will  be  described  in  the  chapter  on  the  computer 
programs.  liquation  (CO),  the  axial  particle  trajectory,  is  identical 
to  equation  (28)  which  describes  the  axial  particle  trajectory  in 
Case  1.  This  is  expected,  because  the  axial  flow  velocity  is 
independent  of  radial  or  tangential  position  and  also  because  the 
force  due  to  the  electrostatic  field  does  not  have  an  axial  component. 
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Case  IV 


The  fourth  case  analyzed  considered  the  motion  of  two  charged 
particles  in  a  free-vortex  flow  field.  Viscous  drag  in  the  fluid 
was  accounted  for  through  the  Stokes  relation,  and  the  particles  were 
assumed  to  influence  the  motion  of  each  other  due  to  the  electrostatic 
body  forces  arising  from  the  electric  charge  on  each  particle.  The 
following  assumptions  were  made  in  this  analysis: 

1)  The  velocity  distribution  of  the  flow  field  is  known; 

2)  The  particle  mass  flow  is  small  compared  to  the  carrier 
fluid  mass  flow; 

3)  The  particles  are  charged; 

A)  The  drag  force  is  given  by  the  Stokes  relation; 

5)  There  is  no  fluid  or  particle  motion  in  the  axial  direction. 

The  equation  of  motion  of  the  first  particle  is  given  by: 

^P,  —/<  =  3tT/<£>  t  t  (61) 

IA 


and  that  of  the  second  particle  is 


pYl^  Ex  >  (62) 

where  and  are  the  electric  fields  at  points  1  and  2  respectively, 
as  shown  in  Figure  7.  Expanding  these  equations  in  (r,9)  coordinates, 
the  following  equations  result: 


V.-,  -  &  -•  isaStfrfr.  +  i  &&  , 

T»i  ff  '  V\ 

Vr.  ^  -t  --  3jta-'/^0  -  V*  )  fr. 


(63) 

(6A) 


-v/r/W  ,  <65> 

<Hrx  n  wi!  '  2 

Vr^—®>  +  7-'- V#1-  --  -  \/0  \  t  ??,  £la  -  (66> 

JA  Vv.  'A  '  1  v 

where  Er^,  Eq^  are  the  radial  and  tangential  components  of  the 
electric  field  at  a  point  (r^,6^)  due  to  the  charged  particle  q^, 
and  E^ »  Eg^  are  t*ie  radifll  a°d  tangential  components  of  the  electric 
field  at  a  point  (r2>®2^  due  to  c^e  charged  particle  q^. 


/fV 


I 


v  / 

V/^-Vr, 

2, 


FIGURE  7 


Referring  of  Figure  7,  Point  has  coordinates  (r^,®^),  a°d 
P_  has  coordinates  (r„,9„). 


The  electric  field  at  due  to  a  charge  at  Point  P  is 


given  by 


Therefore 


Ep'  sine*jy 

ErK  -  x  (coseilL  cose,  t  sin  6, t.  sen e,) 

Ee,  -  ^-4*-  (-coselXsth©, 

£r'-  ^  CDS(e--e'\ 


£e,s  5 tnC©,*-  e>). 

^TT<=a 

The  electric  field  at  P2  due  to  a  charge  qx  at  Point  P]L  is 

given  by  £pS  -^-^(COS  6ia.C  +  5Lh6lXj). 

hTTcoNh  0 


Therefore 


where 


£r*.  a  eiVC0S^r  Stn6,x  cosfr/a-O)  ; 

i  itCo  Hi  ' 

£a  c  (,-C0S6(7.Sth6j.tSChe^Stn(rr/i-G,.'))? 

z  ^r<£o  \ 

Er7  -  41-,  005(0,*- ©xV 

4c5"1(e'’'  e°- 

'K.x^  (r*cos&>.-r;  cos^  +  (r,.sun'3t-r:stne^ , 


"tan 


n,  U09>-  ^  S.tn9i 

r*  cose,.-  n  coso, 


*  tan" 


-i  /  fi.:4ne»L- 

\  f;  OSGx-  b  f  OSG, 
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Therefore  Eri>  Eq  ,  Et^,  and  Eq  can  be  expressed  entirely  in  terms 
of  rx,  Oj,  r2,  and  6 ^ 

Assume  that  the  flow  field  is  an  axisymmetric  free  vortex. 
Therefore , 


Kn  *  o  , 

6(  rz  --  O  , 

(76) 

(X-  04  -  CvAi  ( 

U.0^*  d/ If*,  , 

(77) 

Along  the  particle  trajectories, 

*>■£>  (78) 

Vrt  -  is  ,  Vo,  *  r,  ■  (79) 

At  at 

Substituting  into  Equations  (63)  through  (66),  the  following  equations 
result: 


in  - 

1 

I''  ' 

-  Ci  - 

ISiic 

c*os(e„- 

_e! 

(80) 

cU7- 

Ut, 

I 

dt 

f>vt,4ir6» 

K\ 

r  /<&. 
r,7T^ 

+ 1  — 
At 

it 

w.  <A0.\  _ 

-  71  si) 

Z  sun(e,x-e,s) 

(81) 

^  . 

rJ-jA 

p  in  - 

ii*Je 

C03_(  Ojr- 

(82) 

;U 

1  c«t 

<Wa/lr£a 

a.  - r 

At'- 

+  Z  — 

f  r>. 

Jt 

;  n  ( £ 

rJ(.iv'l3o.-9,') 

- - -x  ,  i  -  — 

7>* 

IX  IT. 

(83) 

where  =  STip-D^/m  and  C2  =  3np.D2/m  . 


I 


These  equations  can  be  nondimensionalized  through  the  use  of 
the  following  parameters; 

X,  --  r,/V0 

'fu*  r>/r“ 

T-  C.-t 

r  £  k  .  r.  C  -  cV, £  -  U-  k  m 

'  c<T*-  \ATj  (AT  qtrcsCi'i.  P 

b(.  il>_  t  ^  K Li  ?7f  7f  •  n,r,'  A.rJ  4a€.gV,  r* 


r0c!Ai  .  ru  K  JACvt-  -  C*G  dj\  + 

ir*  v-VT  /  (AT  4  r  €4  fc>V  a. 


&C,T, il>-T  2  r.  (■,’<£. £!i, gY&  - (87) 
(A'X  AT  AT  Mt?-.  T1,  /  /ic6o^'fl1l,_  i- 


£H‘  -  TTfM't  %  ♦  /1,/k  cc±_(&:S’il  =  o 

d 7 '’  \TT  j  (AT 


X,  f !•  +  d'l'f-J  iV\  t  ?\  -  ^  4  fl,  fi,  *-  ©')  __  0 
6< 1 '  ■  cVT  7  •■•«  /  K<> 

JS.:>  -K  (^A  t  A'.*'-*  -  /l..^  <*>£.(  <>'-<0  -_  0 
C'.T'  'w~  /  c' 1  r 4  '• 


pr 
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p  dQ t  .  dQi/n  di-’i  A  p  \  _  4,  /Is-  Sin  (o*-Qi')  _  n 


(91) 


where  f{x  =  Oc  fo&K*  >  4*  1 


and 


flj-  C,/r.*C,,  /I,*  Cx/C., 

*  '•  /  j±  Stn  9>  ‘  ^  Sl*8>  \ 

6  *  £*-h  V  FUco^e*.-  R>COi0,J 


and  (Kp.  SinGx-  Ri  +  (g*  COS  0a. "  R«  COSQ^ 


Assume  the  following  boundary  conditions: 


at 


d  _  o  £>  ~  r\ 

dT 

d?x  _  <i  Ci  .  r»  d.9x 


T^o,  **  1,  £^-o,  ©,.o,  J&'O, 


^  =  o,  ©*-o,  aa  -o. 

Equations  (88)  through  (91),  subject  to  the  boundary  conditions 
Above,  can  be  solved  using  the  numerical  procedure  mentioned  previously. 
The  computerized  equations  are  given  in  the  chapter  on  the  computer 


programs. 
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NUMERICAL  PROCEDURE  AND  COMPUTER  PROGRAMS 

The  numerical  procedure  used  to  solve  the  equations  governing 
the  particle  trajectories  in  the  radial  and  tangential  directions  is 
the  Runge-Kutta  method  with  Gill's  coefficients.  Initial-value 
problems  of  the  type  encountered  in  this  analysis  can  be  easily 
solved  using  this  method. 

For  a  first  order  differential  equation, 

with  initial  conditions  yQ  and  tQ,  the  basic  formula  in  the  Runge- 
Kutta  method  with  Gill's  coefficients  is  given  by: 

where 

ail,-  *  2«|, 

and 

^  (  J-  "  Ux  +  ( 1  +  ft)  ,  "k*  40 

The  derivation  of  this  method  is  quite  straightforward  and  is 
piesented  in  Reference  49. 

The  following  equations  and  their  associated  boundary  conditions 
were  to  be  solved  using  the  Runge-Kutta-Cill  method. 


■PPPWKKPP 
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Case  I 


££  +  Ji  _  /Cv\’-  6-e~T)T' ,  o 


d?K  „  ai  _  1C2.Y  n-g  >  .  a 
d  T  v~  aT  \W  K' 

ie  _  G  (t-  e‘T)  .  0 

■IT  C,  "R7' 

Boundary  conditions:  T  =  0,  R  =  1,  dR/dT  =  0,  6  =  0. 
Case  II 

(V-  ^iAneJ;  (V-A/r*1,  Cv* 

du  _  Cy  £<  •  Cs-  _L  ,  0 

dr  f  cVf  C,1  R 

Boundary  conditions:  T  =  0,  R  =  1,  dR/dT  =  0. 

Case  III  ,  ,  tNJ-  „ 

HE.  +  -  /<Tf  0~ £~  )  -  ^ 1  -  Q 

.  _  -T  T  IT  l  /-I  1  _  a  /  1  T*> 


,'TV  dT  \c.  y  f?*  r  K 

4  -  (4  <1^3  a  o 

^  t  '  s' 

Boundary  conditions:  T  =  0,  R  =  1,  dR/dT  =  0,  8  =  0. 


It  is  obvious  that  these  three  cases  could  be  solved  by  one 

general  computer  program  through  the  proper  choice  of  the  constants 

C, ,  C„ ,  C. ,  C-  ,  and  Cc. 

1  2  3  h  5 

A  general  program  was  then  written  to  solve  the  following 


equations : 


•T  .{£ i\  (ini  .1  ~  r  c,  •  c-  j_  __  0 


?h  .[ii\ 
dTl  '-T  UJ 


dy  ic 7  \  1-  .  q 

:'T  j  1 

Subject  to  the  boundary  conditions:  T  =  0,  R  =  1,  dR/dT  = 


In  Case  I,  C^  =  C^  =  C,.  =  0.  In  Case  II,  =  0. 
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These  equations  were  reduced  to  a  set  of  first  order  equations 
which  can  be  solved  simultaneously  using  the  Runge-Kutta-Gill  method, 
The  following  transformations  were  used: 

v  F  v, »  «£. 

t<T  IT  ' 

The  three  first  order,  differential  equations  are: 

il=  »  7j 


4.7 


i  y*  v  +.  (£±\ 

«¥  '  "  3  U.j  y,1  '  cP  yk 

cC/<  .  /M  t^elT 
dir  \c*  I  y* 

A  separate  computer  program  was  written  for  Case  IV.  The 
following  equations  and  their  associated  boundary  conditions  were 
solved  using  the  Runge-Kutta-Gill  method. 

£>.lJU  _  TmA^.V  +  dj]  „  A.A>._COS  .  Q 

T 7‘  VJ  i  /  <f  ”  ’ 

Pi  *:G,  +  Ap\  rPtP- e.)  4  y 

V' ■  n  I  ;c,  ■  "'R*’- 


4-,  IL"  -  /!./?«  cv '• 

-,Tr  V“  /  T 


O 


Vv  ^  , 
^  T5' 


.  y  -  +  \  R>\  _  t-LLh  -  :-u>  r 

1  V  ’  .47  y  Ki 


/’4--  /?,•  C'./Cr.’, 

fi«  -  /u,  ^ :  y,  ,  ? 


where 


^  K-j  (- 1  CO  <0|  J 


anu 


,C(06j  i  (:  >.  CGI  9.i  -r\ 


/“  ^  < 
U  vJ  ^ 


The  following  boundary  conditions  apply: 


at 

T-o, 

i/ 

4 

AT 

1 

0  i  :  0; 

d&, 

'  dJ 

If 

O 

■Rs.: 

■ 

d  .9  . 
AT 

-  o. 

0a=  O, 

_ 

d  T 

G. 

The  following  transformations  were  used  to  obtain  the  first  order. 


differential  equations: 


Y*  - 


Y®*- 

A  Vs. 

AT 

AT 

Y«  = 

<?, 

Yr  : 

An, 

A  Y« 

AT 

AT 

Y  = 

Yl  ' 

A  •  i 

A  Vo 

AT  ’ 

'  AT 

Yfc  •  0s. 


Y, 


dT 


dVv 

‘ciT 


The  first  order  equations  which  were  then  programmed  are: 


=  y. 

c'T 


Ta  i  y.  ’  -  Yj  -  ,/S,  cci(e-  y.  W 

at 


a  V*  _  y. 


(  T 


r/s-  -  - >V ('  j? V?  //  +  r)  +  4j/V>  -  /i>  V0  C<  n  (g-VY/Y/,. 


Cl 


where 


11’ 

dT 

P  z  X/V^-  ^/Vt  +  COS(e-y«,)/R‘t' 

y, 

A>_.  _  _  Vq  (7V-,  /vc  «-  r  /Ij^-v/Vb^4  J- C-i  n  (0  -  'lO/V  Yt 
dT 

-i  /  yt  sm  v t  -  v>  cc^v*. 

<0-  tan  (  y 0  :div4  -  Vs  cosy*, 


-  V, 


and 


k^z-  (/t  sen  Vt  -Y> SimV-j')  +  (Vt  coi>  Vi  -  V?  y«)  . 


The  two  computer  programs  used  to  determine  the  particle 
trajectories  for  each  of  the  five  cases  considered  are  listed  here. 
The  computer  program  PARVOR  was  used  in  the  determination  of 
particle  trajectories  for  the  first  three  cases.  The  computer 
program  TWOPAR  was  used  to  analyze  the  motion  of  the  two  charged 
particles  in  Case  IV. 
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Y'f-  1  Wu  J  (I  ^ill^'Llrull 

I’  ►  Nl'  .  >  '  /  i  »  i  t  /  (  (Art  I  .U  (,-  I  )  .(  ,■  uv  I  ,C'.  {/>.»»  .(,'■>  <  in)  .  A  ( /> 

-I  *■  A  *  I  It  1,  s'  t  .  ••  t  S 

1  f  1  (  l  «  ,  /  ||  ,  ) 

M  *  ft  I  t  .  t  -  i  • • 

Mlti  1  /*«{,*,  ►  I  i  |  ”  M  L 

?  f  *' v  .  ..  i  1 1  ».?..!  .  i 

r*  *  •’«  |  "  P  t  1  ll')t-V=l*M)) 

-ft  ft’>  r;.IUI'.|tr:h'/l 
I'*  t.  j  .  (■  H'  :  I  r  tl 

tV  ►  -til  P  t  1 1.  -  ( **  I  t  f,  =  |  t  f  j  S  1 

r.'  =  *<  •  1 

i ,  ^ f  1  s  ,  ♦  i 

ri  ft  p.  it  1 1  )  ,  |- 1  .;.m 

i  n  i  j  |  =  1  .*i* 

n  * ( I ) =r ( I ) 

l"-  II  1 1  =  1  .'.I 
|l'»  1  I  L^-  I  t'1’ 

1)1  >  II  L  '•  s  I  « ' 1 '« 
lj.>  I  I  l»  ft  =  I  • '  " 

c^ssO'd'’  i  /n  a  1 1 

C<1  *  I  sCJ’tf.  II  1  “CM  (|.<>)  /c  1  (L  I  1  •*> 

.  i  jti,tM'Ht(,iiLl),C£lL'’itCJiM.ii.  *•>.':=.  (i  ft)  .C0rj.CoM  *»kuT  .  f“<«< 
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CKOM  I  =  ».t  I  1.1/  I  <>%  .*l>KL  t  =  «  «K  1  i.b.J-x  .  *»  |Mft*s<.  ,F.  H.m 
r>r  M'j  I  = 

1  '■  ft  =  . 

Y  ■>  (  I  )  =  I  . 

Y  ft  (  P )  r  i  (  .1 ) 

TP(  llni.,'  <-  ?»  (I.  -  r»P<-Y(i)))<''-V/Ylc'»'»**J.Cii\l/Yl?>-Yt  '• 

Yi>(  )  =1  i  .<■  (  1  . <o  (  -  Y  (  I  )  J  )  /  Y  <C  )  «“? 

AVffl.  =  r  (  S)  «  i  ■>  (ft  | 

.V  n5i  (  •  I  “  i  >i  ./  i.  1  •»  |  Mm 

I'PIM  . 

4  )  *  H  1/ I  M*  .  ■>  I  |'->  <>  .  ■  •  «hUlllOS*'  »*><»»MAl)l  ftl.  Vf.L1*  «?<♦<* ‘H-OLt  (k*»0)«*<?a 
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P  *'  1  ’ :  l  it  I  Y  I  I  I  •  1  =  1  .  '  ' .  )  .  A  V  .  L  ♦  ft  M  o 
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(i  1  =  1  t  I 

h  ft’-c  (  1  )  =  •. . 
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A--C  (  '  •  ft  I  =  . ft';ft,|M?ft  •  I.H  1  <t 

ftMl  |  •  l|  :l  .  I  'll  |»  /Ml  lfr*Y 

ftl'C  I  '  ■  “»  1  =  •  I  1 1 11  h  t>  -  >  ll  I;  M  fc  / 

A » •  (  1  ;i* 

ft*  C  l  •'  “’1  =  1. 

/••'  C  (  ♦/>  =  !. 
tf'C  ('  :•  I  =  •  . 
ft  i>C  (  *. 

/••  C  (  •  I*  1  )  =  .  ft'A  /.V  I  »'*i  I  ft 
ft- iC  <  •  ♦  1  i  >  =  I  -  /  /t'J«. /f'ltrtft 

ftnC  (  •!<:)  =  .' 

7  |  i'  o  j  =  .  i , : i'  i 

I  .  ft  *  '1  1  =  1  til', 

=  1  J*  I  >  -  (  l  I  I  -  Al-C  U*b)  I’A-lCI  I  )  ) 


l>-  .  .  Tv',it*.,  n  |  >  .t-i  1 1  ,  ')t_  1  n  1 1  ) 

•  x  ! •"  i  (  : )  •  i • 1  t  i)  i>i  i^,!l 

1  • » • 

1  b  i  •  T  (  I  .  !\  '  I 

►  •  >  2  •  i  •  • . .  •  li  r  •  ('  L  l  1  »  t  '■ "  < 

?  h  •  <  •  til  .  .  /•  I  .  - ) 

P '  1  1  I  <i.l  •  '  '  • « .  a  <• .  a  S  ,  i , :  l  |  <  l  i.  < 

I  b  1  I  1  !  f'  >  •  ’>  I  s1*  i  i  ►  »  I  ^  »  •“Alsfc  »^  |  (*  )  t  (n  i  o.'i4s((» 

C  r.  I  *  •  *  »(*  \  ^  =  w  •  *  !  *  ,  s  /  1 2  *  i  <*  i  j  c.  C.  *  z  '*  «  K  &  %  °  7  A  A  s  •  1. 1  §  0  ) 

V  =  > 

:.r  \s  1  a 

mi  ;.i>  /  •  ( t  c  i  >  ,  i  s  i > 

Kl  I  II 
1  1  A  =  . 

TP  (  I  )  =  I  . 

Y  b  (  X  )  =  <  (  J  1 

>ri)  r:(|i  i»r  /,  ( i  ( r>) )  -  y  ( 2 1  «c  (j  ( >  ( 4 1  ) 

rCn r  (’• )  <>M  t  r  t«>)  )  -  »  (  2)  «>s  |M  r  ( )  | 

!»t  i,t..|  (  ft;  i  .  /  ,  n,  ) 

«S  t  s  •  C  1  *  1  ■  *»  a  f*  »■  «  A  L  Y  C 

YP  7 .71  =  T  I  r  1  vy  (  ,1  -  Y  (  S)  -  Y  (  j)  -A  1  “  (  IMt.T  A-Y  (<.))/(MbT*-<ST> 

YP  ( )  =  r  (  : ) 

yP(i)s-y(M*"(p.*'Y('})/Y(2)*I.)»aJ/(y(?i*v(/'))-a1«a2*S'  y(T.<Ma- 

-  Y  (  A  )  I  /  I .  M  « . '  ,  1  <•  i  (  2  )  ) 

YP(/)  =  f(<i)«Y(  >)  «*  Y  (  J)  -AhI'Y  »  n  1  7M(.  |  A-Y  C  H)  )  /  it  at*' MSI) 

( ■. ) 

Yp  1 1. )  r._  T  I  ;  I  •'  (  2.  **Y  (  I)  /  Y  (  Si  «  «*«  )  ♦  A  |n.i4/|YlM*»tM  )  ♦  \  1«A  ,“S  1  -  I  I mE7  A- 

-  Y  (  *  I  )  /  (  -  S  1  *'  ^  I  •'  V  (  fi )  ) 
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t  ( ;  i  -  y  i  i  i  •  ■  !  t  i  ” ' r . 
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M  Af  C  (  l  )  =<*"C  l  I  )  »  i.  «  I  ►  -0-AHC  (.)♦>#»  "YP  t  |  ) 

#*•(11*1. 
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RESULTS  AND  DISCUSSION  OF  RESULTS 

The  results  presented  and  discussed  here  are  based  on  the  data 
presented  in  Appendix  B.  The  numerical  procedure  and  computer 
programs  discussed  in  the  previous  section  were  used  to  calculate 
these  results. 

For  the  case  of  an  uncharged  particle  moving  in  a  free-vortex 
flow  field,  Case  I,  particle  trajectories  have  been  determined  for 
various  values  of  particle  density  and  diameter,  and  various  values 
of  free-vortex  flow  velocity.  Figures  8  through  13  describe  the 
particle  trajectories  in  the  (r/rQ,0)  plane  and  the  radial  migration 
as  a  function  of  dimensionless  time.  The  axial  particle  trajectory 
as  a  function  of  dimensionless  time  is  plotted  in  Figure  14. 

Figures  8  through  11  indicate  that  as  the  particle  density  and 
diameter  increase  at  constant  free-vortex  velocity,  the  entrainment 
of  the  particle  in  the  free-vortex  is  decreased.  Figures  12  and  13 
indicate  the  radial  migration  of  the  particle  tends  to  increase  as 
the  free-vortex  flow  velocity  increases  while  the  particle  density 
and  diameter  are  held  constant. 

These  results  can  be  explained  by  examining  the  parameter 
2 

(Cv/CrQ  ),  which  represents  the  ratio  of  the  centrifugal  or  inertial 
force  per  unit  mass  on  a  particle  moving  with  a  velocity  equal  to 
the  tangential  velocity  of  the  flow  field  to  the  drag  force  per  unit 
mass  on  the  particle.  As  such,  the  parameter  is  indicative  of  the 
relati\e  importance  of  the  centrifugal  or  inertial  force  compared 
with  the  drag  force  for  a  given  particle.  As  the  particle  density 


r 


or  diameter  increases,  the  parameter  (Cv/CrQ  )  increases,  indicating 
that  the  centrifugal  force  tends  to  dominate,  resulting  in  an 
increase  in  the  radial  migration  and  conversely,  reducing  the 
particle  entrainment.  As  the  free  vortex  flow  velocity,  represented 
by  the  parameter  Cv ,  increases,  the  centrifugal  force  once  again 
tends  to  dominate,  increasing  the  radial  migration. 

The  results  also  indicate  that  the  angular  displacement  of  the 
particle  is  strongly  dependent  on  the  free-vortex  velocity.  Examination 
of  Figure  13  indicates  that  as  the  free-vortex  flow  velocity  is  in¬ 
creased,  radial  particle  migration  increases,  and  the  angular  dis¬ 
placement  of  the  particle  increases  also.  However,  at  constant 
free-vortex  velocity,  as  the  radial  particle  migration  increases,  the 
angular  displacement  of  the  particle  decreases.  Figures  9  and  11 
illustrate  this  behavior.  These  results  can  be  explained  by  examining 
the  differential  equation  governing  the  angular  displacement  in  Case  I: 


d«  _  „  ..  -Ctw  2 

Tt  cv(1  '  e  )/r 


The  angular  displacement  is  seen  to  be  directly  proportional  to  the 
free  vortex  flow  velocity  parameter,  C^,  and  inversely  proportional 
to  the  square  of  the  radial  position  of  the  particle. 

Examination  of  the  equation  governing  particle  motion  in  the 
axial  direction  indicates  that  it  is  completely  decoupled  from  the 
equations  of  motion  in  the  radial  and  tangential  directions.  This 
will  be  true  if  the  axial  flow  field  is  assumed  to  be  independent  of 
the  radial  and  langcntinl  coordinates  of  the  particle.  This  is,  of 
course,  the  assumption  made  in  these  analyses.  Figure  14  indicates 
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that  as  the  axial  tluid  velocity  becomes  more  negative,  the  particle 
is  displaced  a  larger  axial  distance,  and  the  terminal  velocity  of 
the  particle  becomes  larger,  as  expected.  The  axial  fluid  velocity 
in  effect  controls  tno  settling  rate  of  the  particle,  and  as  such 
would  be  an  important  parameter  in  a  centrifugal  separator,  since  it 
would  be  desirable  to  obtain  separation  of  particles  prior  to  their 
settling  out  of  the  system. 

The  second  case  considered  was  that  of  a  positively  charged 
particle  moving  in  a  viscous  medium  under  the  influence  of  an  applied 
electric  field  due  to  a  positively  charged  line  source.  Radial 
particle  trajectories  as  a  function  of  dimensionless  time  are  shown  in 
figures  15  through  18  for  various  values  of  particle  density  and 
diameter,  line  source  strength,  and  particle  charge- to-mass  ratio. 
Figures  15  and  16  indicate  that  as  the  particle  density  and  diameter 
increase  at  constant  charge- to-mass  ratio  and  line  source  strength, 
the  radial  migration  increases.  Figures  17  and  18  indicate  that  as 
the  particle  charge- to-mass  ratio  increases,  or  the  line  source 
strength  increases  while  particle  density  and  diameter  are  held 
constant,  radial  particle  migration  also  increases.  As  the  particle 
density  and  diameter  increase,  the  drag  force  per  unit  mass  is 
decreased,  permitting  an  increase  in  the  radial  particle  migration. 

As  the  line  source  strength  or  particle  charge- to-mass  ratio  is 
increased,  the  repulsive  force  between  the  charged  particle  and  t lie 
line  souree  is  increased,  since  both  arc  assumed  to  bo  positively 
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charged.  Tlic  result  is  an  increase  in  the  radial  particle  migration 
also . 

The  third  case  analyzed  the  trajectories  of  a  positively  charged 
particle  in  a  free-vortex  flow  field,  subject  to  an  applied  electric 
field  due  to  a  positively  charged  line  source.  This  case  in  essence 
combined  the  forces  acting  on  the  particle  which  were  considered 
in  the  two  previous  cases.  Figures  19  through  28  illustrate  the 
particle  trajectories  and  radial  particle  migration  as  a  function  of 
dimensionless  time  for  various  values  of  particle  density  and  diameter, 
free-vortex  flow  velocity,  line  source  strength,  and  particle  charge- 
to-mass  ratio.  The  trends  in  radial  particle  migration  as  a  function 
of  these  parameters  which  were  observed  in  Casps  1  and  11  were  ob¬ 
served  here  also. 

In  order  to  assess  the  relative  effects  of  each  of  the  forces 
involved,  centrifugal,  drag,  and  electrostatic,  this  analysis  was 
repeated  assuming  the  particle  was  negatively  charged.  This  assump¬ 
tion  results  in  the  electrostatic  force  and  the  centrifugal  force 
opposing  each  other,  rather  than  reinforcing  eacli  other  as  in  Case  III. 
Figures  29  through  38  give  the  particle  trajectories  and  the  radial 
particle  migration  as  a  function  of  dimensionless  time  for  Case  11]  A. 
Figures  29  tl  rough  32  indicate  that  at  constant  particle  charge— lo- 

raass  ratio,  line  source  strength,  and  free-vortex  velocity,  as  tin 

* 

particle  density  and  diameter  increase,  radial  migration  reaches  a 
limiting  radius  at  which  the  particle  becomes  entrained  in  the  free- 
vortex  The  attractive  electrostatic  force  between  the  positively 
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charged  line  source  and  the  negatively  charged  particle  and  the 
drag  force  in  effect  tend  to  balance  the  centrifugal  or  inertial 
force,  allowing  the  particle  to  reach  an  equilibrium  position. 

Figures  33  and  34  indicate  that  as  the  free-vortex  flow  velocity  is 
increased  with  all  other  parameters  held  constant,  radial  particle 
migration  is  increased.  This  indicates  that  the  centrifugal  force 
becomes  dominant  as  the  free-vortex  flow  velocity  is  increased. 
Examination  of  Figures  35  and  36  indicates  that  as  the  line  source 
strength  is  increased  with  all  other  parameters  held  constant, 
particle  entrainment  is  increased.  Thus  the  attractive  electrostatic 
force  between  the  positively  charged  line  source  and  the  negatively 
charged  particle  tends  to  balance  the  centrifugal  force  on  the 
particle.  This  same  trend  is  noted  as  the  particle  charge  to  mass 
ratio  becomes  more  negative,  as  evidenced  by  the  results  shown  in 
Figures  37  and  38. 

Examination  of  the  radial  particle  trajectories  shown  in 
Figure  29  indicates  that  at  times  shortly  after  the  particle  is  re¬ 
leased,  its  motion  is  inward  toward  the  center  of  the  free-vortex. 
This  can  be  explained  as  follows:  since  initially  the  particle  is 
motionless,  some  time  is  required  for  it  to  achieve  a  tangential 
velocity  of  sufficient  magnitude  such  that  the  centrifugal  force  is 
larger  than  the  electrostatic  force.  Until  this  occurs,  the  net 
motion  of  the  particle  will  be  inward.  It  is  also  evident  that  the 
particle  having  the  greatest  inward  penetration  will  achieve  an 
equilibrium  radius  earliest. 
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In  order  to  assess  the  relative  effects  of  the  forces  considered 

to  act  on  a  particle  released  in  a  free— vortex,  comparisons  between 

the  three  cases  just  discussed  were  made  and  the  results  presented  in 

Figures  39  through  44.  Figures  39  and  40  compare  the  radial  particle 

migration  and  particle  trajectories  for  Case  1,  an  uncharged  particle 

entrained  in  a  free-vortex,  and  Case  11,  a  charged  particle  moving  in 

a  viscous  medium  under  the  influence  of  an  applied  electric  field. 

It  is  evident  that  the  repulsive  electrostatic  force  lias  a  greater 

effect  on  radial  migration  than  the  centrifugal  or  inertial  force. 

This  is  due  to  t lie  fact  that  the  electrostatic  force  is  proportional 

3 

to  1/r),  whereas  the  centrifugal  force  is  proportional  to  (1/r  ). 

Thus  greater  radial  migration  is  obtained  from  a  particle  moving  under 
the  influence  of  the  applied  electric  field.  In  Figure  40,  the 
particle  trajectories  for  these  two  cases  a -e  compared.  As  expected, 
the  particle  trajectory  for  Cast  II  is  a  straight  line,  since  in  this 
case  the  particle  has  no  initial  tangential  component  and  the  only 
forces  acting  on  it,  electrostatic  and  drag,  have  radial  components 
only . 

Figures  41  and  42  provide  a  comparison  of  the  results  from  Cases 
II  and  III.  It  is  seen  that  inclusion  of  the  centrifugal  force  due  to 
the  presence  of  the  f roc-vort c x  results  in  an  increase  in  radial 
migration  in  Case  Ill  over  that  predicted  for  Case  1J.  The  centri¬ 
fugal  force  augments  the  radial  electrostatic  force  ,  increasing  radial 
part  ic  le  migration. 
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Figures  43  and  44  compart-  the  results  from  Case  1,  Case  111,  and 
Case  III  A.  When  the  line  source  and  the  particle  are  both  positively 
charged,  as  in  Case  111,  the  electrostatic  force,  which  is  repulsive 
in  nature,  augments  the  centrifugal  force  and  causes  the  radial 
particle  migration  to  increase.  When  the  line  source  and  the  particle 
are  assumed  to  he  oppositely  charged,  as  in  Case  I'll  A,  the  results 
indicate  that  the  electrostatic  force,  which  is  now  attractive,  and 
the  centrifugal  force  tend  to  balance  each  other  out,  allowing  the 
particle  to  attain  an  equilibrium  radius.  The  particle  trajectories 
shown  in  Figure  44  also  exhibit  the  same  behavior  noted  previously  -- 
as  the  radial  migration  increases,  the  angular  distance  the  particle- 
travels  is  decreased. 

The  fourtli  case  considered  analyzed  the  motion  of  two  charged 
particles  in  a  free— vortex  flow  field  under  the  influence  of  the 
electric  field  due  to  the  charge  on  each  particle.  Calculations  were 
made  to  determine  the  effect  on  particle  motion  of:  particle  charge- 
and  the  sign  of  the  charge,  and  particle  mass.  The  particle  tra¬ 
jectories  and  radial  particle  migration  are  presented  in  Figures  45 
and  44.  The  results  indicate  that  field  effects  due  to  particle 
charge  on  particle  motion  arc  negligible,  if  the  particle  charge  is 
of  the  same  magnitude  as  that  considered  in  Cases  |]  and  ill. 

The  results  also  indicated  that  changing  the  sign  of  Ihe  charge  on 
one  of  the  particles,  in  if  feet  changing  the-  electrostatic  force  from 
a  repulsive  force  to  an  at  tract  ivi  force,  had  no  eliccl  on  particle 
motion.  Increasing  tin  mass  of  the  particle  assumed  to  It  initially 
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farther  from  the  center  of  the  free  vortex  increased  its  radial 
migration,  since,  as  expected  from  ttie  results  obtained  in  Cases  1 
and  Ill,  the  drag  force  acting  on  the  particle  was  reduced.  Since 
the  parameters  governing  the  motion  of  the  particle  nearer  the  center  of 
the  free  vortex  were  not  changed,  its  motion  was  unchanged  from  the 
previous  prediction.  An  increase  in  the  magnitude  of  the  particle 
charge  by  a  factor  of  2000  did  result  in  a  substantial  change  in  the 
radial  particle  migration  and  the  particle  trajectories.  The  repulsive 
electrostatic  force  between  the  charged  particles  caused  the  particles 
to  separate  rapidly  before  the  frcc-vortex  flow  field  began  to  aifcct 
their  motion.  Relatively  large  particle  charge- to-mass  ratios  are 
required  before  field  effects  due  to  the  charge  on  the  particles 
would  he  of  sufficient  strength  to  affect  particle  motion,  A 
particle-flow  system  containing  a  large  number  of  charged  particles 
may  result  in  electrostatic  fields  of  sufficient  strength  to  affect 
the  particle  motion 

Any  extension  ni  this  analysis  should  attempt  to  account  for 
field  effects  on  fluid  motion.  In  the  one  and  two  particle  systems 
considered  here,  field  effects  on  fluid  motion  are  no  doubt  small 
due  to  the  highly  localized  body  force  acting  on  the  fluid  as  a 
result  of  the  electrostatic  field.  This  may  not  he  the  ease  in  a 
particle  system  composed  of  a  large  number  of  charged  particles,  in 
which  the  electrostatic  body  force  may  cause  the  fluid  to  be 
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CONCLUSIONS 


lor  the  case  of  an  uncharged  particle  moving  in  a  frec-vortex 
system,  the  results  indicate  tliat  at  constant  frec-vortex  velocity, 
entrainment  of  a  particle  is  decreased  or  conversely,  the  radial 
particle  migration  is  increased.  Radial  particle  migration  also 
is  increased  as  the  frec-vortex  flow  velocity  is  increased  while  the 
particle  density  and  diameter  are  held  constant.  As  radial  particle 
migration  increases,  the  angular  displacement  of  the  particle  is 
decreased.  Since  the  axial  flow  velocity  was  assumed  to  be  constant, 
the  equation  of  motion  of  the  particle  in  the  axial  direction  is 
independent  of  the  radial  and  tangential  coordinates  of  the  particle. 
Thus  particle  motion  in  the  axial  direction  is  not  influenced  by 
what  is  occurring  in  the  radial  and  tangential  directions. 

For  the  case  of  a  charged  particle  moving  in  a  viscous  medium 
under  the  influence  of  a  charged  line  source,  particle  motion  is 
strongly  influenced  by  the  particle  charge- to-mass  ratio  and  the  line 
source  strength. 

For  the  case  of  a  positively  charged  particle  moving  in  a  free 


vortex  under  the  influence  of  an  applied  electric  field  due  to  a  line 
source  of  po.sitivc  charges,  the  electrostatic  force  reinforces  the 
centrifugal  force,  causing  an  increase  in  radial  migration  over  the 
cases  In  which  each  force  is  considered  to  aet  individually.  If  the 
particle  is  assumed  Lo  be  negatively  charged,  the  electrostatic 
force  and  the  centrifugal  force  tend  to  balance  each  other  out, 
allcuvin  the  particle  to  attain  an  equilibrium  radius.  The  drag 
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force  alone  is  not  of  sufficient  magnitude  to  balance  either  the 
electrostatic  force  or  the  centrifugal  force. 

Results  from  the  case  of  two  charged  particles  moving  in  a 
frec-vortcx  indicated  that  the  electric  fields  due  to  each  of  the 
particles  were  not  strong  enough  to  affect  the  motion  of  the 
particles.  If  a  system  with  a  large  number  of  charged  particles  were 
to  be  analyzed  in  a  similar  manner,  electric  field  effects  due  to 
the  particle  charges  would  influence  the  particle  motion. 
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RECOMMENDATIONS 

'flic  following  recommendations  arc  made  as  a  result  of  this 
study : 

1)  The  analysis  should  be  repeated  using  a  more  complex 
vortex  system  with  a  more  realistic  axial  flow  distribution.  Some 
of  the  more  complex  vortex  systems  mentioned  in  the  literature 
survey  could  be  incorporated  in  this  analysis. 

2)  The  analysis  should  be  extended  to  a  system  incorporating 

a  large  number  of  particles,  with  the  particles  preferably  distributed 
in  size,  so  that  the  vortex  could  be  used  as  a  particle  separator. 

A  charge  distribution  should  also  be  imposed  on  the  particles, 
so  that  electrostatic  effects  could  be  studied. 

3)  Field  effects  on  fluid  motion  should  be  determined  to  see 
if  this  would  have  a  significant  effect  on  particle  motion.  It  is 
believed  that  in  a  one  or  two  particle  system,  field  effects  on  fluid 
motion  are  negligible. 
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Appendix  A 


nomenclature 


Svnil  o  1 
B 


C 


C 


D 


C 


V 


D 

P 

E 


F 

g 


A 

i 


j 


m 

P 


r 


R 


t 

T 

u 

V 


z 


Z 


P 


_ Descript ion _ 

Body  force  net  in;;  on  the  particle 

Particle  drag  parameter 

Particle  drag  coefficient 

Free-vortex  constant 

Particle  diameter 

Electrostatic  field 

Total  force  actin;;  on  tilt'  particle 

Acceleration  due  to  gravity 

Gravitational  constant 

Unit  vector  in  the  x-direction 

Unit  vector  in  the  y-direction 

Particle  matis 

Particle  charge 

Radial  coordinate 

Dimensionless  radial  coordinate,  R=r/r 

o 

Time 


_ Un  i  t  s 

lbf 

-1 

sec 

2 

ft  /sec 
ft. 

lb^/coul 

lbf 

32.2  ft/ sec 2 

32.2  fl-lb  / lb. 

m  i 


lb 

m 

coul 

ft 

sec 


Dimension]  i  as  time,  T-CL 

Flow  field  velocity  ft/sec 

Particle  velocity  ft /sec 

Axial  coordinate  ft 

Dimensionless  axial  coordinate,  V~-v.lv. 

o 

PI  u  i  d  le  ns  i  i  y  ]  h  /  f  L  ^ 

m 


2 

-  sec 
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Svmbo  1 

Description 

Units 

Particle  density 

lb  /ft3 
m 

Fluid  viscosity 

lb  /ft- sec 
m 

8 

Tangential  coordinate 

- 

X 

Line  source  strength 

coul/ ft 

c 

o 

Permittivity  constant 

„  ,,  ,,.-12  coul 
3.66x10  _ 

1L.  ^ 

lb{-fi 


Subscripts 

D 

P 

r 

G 


Description 

l>rag 

Particle 

Radial  component 
Tangential  component 


z 


Axial  component 
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Appendix  B 

DATA 

The  data  used  in  the  analysis  of  each  of  the  four  cases  which 
have  been  discussed  will  be  presented  here.  To  determine  the 
effects  on  particle  motion  of  variations  in  particle  density  and 
diameter,  and  free-vortox  flow  velocity,  these  parameters  were  varied 
through  the  following  range  of  values: 

Particle  diameter,  D  ,  50,  100,  200  microns 

3 

Particle  density,  p  ,50,  200  lb  /it 
}  p  m 

2 

Free-vortex  flow  parameter,  C ,  20,  50,  200  ft  / sec 

To  determine  the  effects  on  particle  motion  of  charged  line  source 
strength  and  particle  charge- to-nass  ratio,  these  parameters  were 
varied  :in  this  analysis  through  the  following  ran"e  of  values: 

Line  source  strength,  X,  1.03x10  1.03x10  \  1.03x10  ^  coul/ft 

Particle  charge- to- mass  ratio  (q/m)  4.54x10  4.54x10  , 

4 . 54x10'**  coul/lb 

in 


